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SUMMARY 


Under  the  assumption  of  a  Stationary  Parameter  Long  Observation  Time 
(SPLOT)  process,  the  optimum  multisensor,  multitarget  time  delay  processor 
was  derived  from  a  Maximum  Likelihood  (ML)  viewpoint.  The  resulting  processor 
was  obtained  by  reducing  the  vector  likelihood  equation  via  straightforward, 
but  somewhat  tedious,  manipulations  to  the  simplest  form.  The  optimum  multi- 
sensor  multitarget  processor  provided  the  basis  for  the  subsequent  detailed 
studies  on  time  delay  vector  estimation,  localization  parameter  estimation, 
power  spectral  estimation,  suboptimum  processor  realization  and  post 
Generalized  Cross  Correlation  (GCC)  multitarget  parameter  estimation.  Major 
findings  of  this  study  are  summarized  below. 

1.  The  optimum  multisensor,  multitarget  time  delay  estimator  is  a 

highly  coupled,  multichannel  signal  processor.  The  estimator's  per¬ 
formance  in  terms  of  the  Cramer-Rao  Lower  Bound  (CRLB)  was  evaluated 
for  the  following  cases:  one  target  and  two  sensors,  two  targets  and 
two  sensors,  and  one  target  and  M  sensors.  For  the  one  target,  two 
sensor  case,  the  optimum  multisensor  multitarget  processor  reduces 
to  the  GCC  studied  by  Knapp  and  Carter  (Reference  1).  For  the  two 
sensor,  two  target  case,  the  resulting  processor  was  studied  in 
detail.  A  closed  form  analytical  expression  was  obtained.  The 
results  of  the  study  showed  that  the  optimum  processor  provides  a 
significant  improvement  over  a  conventional  GCC  processor.  For  the 
one  target,  M  sensor  case,  the  matrix  CRLB  was  obtained. 

The  results  of  this  study  indicated  that  for  an  M  sensor  array,  the 
M-l  inter-sensor  time  delays  can  be  obtained  from  M-l  correlations. 
This  is  a  significant  reduction  from  Hahn's  approach  (Reference  25), 
where  a  total  of  M  (M-l)/2  correlations  are  required.  Additionally, 
the  results  showed  that  the  variance  of  the  time  delay  estimate 
between  any  two  sensors  decreases  with  M,  the  total  number  of 


sensors. 


The  localization  parameter  estimation  was  examined  in  detail  for  a 
three-sensor  array  from  two  perspectives:  (a)  a  one-step  focused 
beamformer  approach  where  a  direct  range  and  bearing  estimate  is 
sought,  and  (b)  a  two-step  time  delay  to  range  and  bearing  approach 
where  time  delay  estimates  are  first  sought  and  range  and  bearing 
estimates  are  obtained  via  a  geometric  mapping.  The  study  indicated 
that  both  approaches  yield  an  identical  performance  bound.  However, 
for  practical  considerations,  the  two-step  approach  is  generally 
preferred  because  of  the  symmetrical  property  of  the  GCC  function  in 
the  time  delay  variable.  Furthermore,  the  localization  performance 
based  on  the  optimum  time  delay  processor  was  compared  to  the  conven¬ 
tional  approach.  The  result  of  this  study  showed  that  the  optimum 
processor  yields  a  one-sigma  localization  error  ellipse  which  is 
approximately  one-half  smaller  than  the  conventional  approach.  This 
improvement  comes  directly  from  a  better  bearing  estimation.  The 
range  variance  is  identical  between  the  optimum  and  the  conventional 
approach. 

A  major  assumption  used  in  the  derivation  of  the  optimum  multisen¬ 
sor,  multitarget  time  delay  processor  is  the  known  target  and  noise 
spectra.  To  relax  this  constraint,  a  joint  time  delay  and  spectral 
estimator  is  derived.  The  result  of  this  study  showed  that  the  time 
delay  estimate  and  the  spectral  estimate  are  uncorrelated.  This 
implies  that  for  the  unknown  target  spectrum  case,  a  joint  time  delay 
and  spectral  estimator  can  be  implemented.  The  spectral  estimation 
process  does  not  degrade  the  performance  of  the  time  delay  estimate. 
Furthermore,  it  was  found  that  while  the  variance  of  the  time  delay 
estimate  decreases  as  an  inverse  function  of  the  observation  time, 
the  variance  of  the  spectral  estimate  decreases  as  an  inverse  square 
function  of  the  observation  time. 

The  optimum  multisensor,  multitarget  time  delay  processor  is  an 
order  of  magnitude  more  complex  than  the  conventional  GCC  processor. 
Therefore,  for  a  practical  implementation,  a  suboptimum  realization 
should  be  considered.  One  suboptimum  procedure  is  to  assume  a  low 
target  signal-to-background  noise  environment.  It  was  found  that 
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the  resulting  processor  is  significantly  simplified.  From  the  mul¬ 
titarget  viewpoint,  the  conventional  GCC  processor  can  be  considered 
as  a  suboptimum  processor.  The  performance  of  the  GCC  processor  is 
compared  to  the  optimum  processor  as  a  function  of  time  delay 
separation. 

5.  Because  of  its  simplicity,  the  conventional  GCC  processor  is  imple¬ 
mented  in  many  existing  sonar  systems.  On  the  other  hand,  the  appli¬ 
cation  of  the  optimum  multisensor,  multitarget  processor  requires  a 
major  modification  of  the  conventional  GCC  approach.  Therefore,  an 
alternate  approach  is  to  provide  additional  multitarget  processing 
at  the  GCC  output.  For  this  purpose,  a  post  GCC  multitarget  estima¬ 
tor  (or  more  appropriately  the  matched  estimator)  was  investigated. 

In  essence,  the  matched  estimator  determines  the  best  estimate  of 
the  unknown  parameter  vector  from  a  reference  function  which  matches 
the  observed  noisy  GCC  output  under  a  Least  Mean  Square  (LMS)  cri¬ 
teria.  The  matched  estimator  was  simulated.  The  simulation  results 
were  compared  to  the  theoretical  predictions  as  well  as  to  the  opti¬ 
mum  processor.  The  results  of  this  study  indicated  that  the  matched 
estimator  provides  comparable  performance  with  the  optimum 
processor. 
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OPTIMUM  MULTISENSOR,  MULTITARGET  TIME  DELAY  ESTIMATION 


CHAPTER  1 
INTRODUCTION 


1.1  BACKGROUND 

In  passive  sonar  signal  processing,  one  of  the  basic  measured  variables 
is  the  time  difference  of  arrival  (TDOA),  or  simply  the  time  delay  between  two 
sensor  arrays.  Since  the  time  delay  vector  from  targets  to  sensors  is  func¬ 
tionally  related  to  the  target-sensor  geometry,  the  final  parameters  of 
interest  (i.e.,  the  target  location  parameters)  can  therefore  be  obtained 
(for  appropriate  sensor  distribution)  via  a  direct  inverse  relation  from 
measured  time  delays. 

The  measurement  of  time  delay  is  mechanized  through  a  Generalized  Cross 
Correlation  (GCC)  function.  The  GCC  is  derived  as  a  Maximum  Likelihood  Esti¬ 
mator  (MLE)  operating  between  any  two  sensor  arrays.  Detailed  discussion  on 
the  subject  of  the  GCC  can  be  found  in  Knapp  and  Carter  (Reference  1),  Carter 
(Reference  2),  Hahn  and  Tretter  (Reference  3),  and  Hassab  and  Boucher  (Refer¬ 
ence  4'. 

In  the  above  literature,  the  GCC  is  optimized  under  a  single  target 
assumption  (or  more  specifically,  a  single  coherent  noise  source).  In  the 
presence  of  multiple  targets,  or  multipath  environment,  there  are  multiple 
correlation  peaks.  The  existence  of  multiple  correlation  peaks  causes 
performance  degradation  to  the  existing  measurement  system.  The  extent  of 
this  degradation  is  a  function  of  signal  spectral  characteristics,  signal-to- 
noise  ratio  (SNR),  signal-to-interference  ratio  (SIR),  and  the  relative  time 
delay  separation. 

In  a  multisensor,  multitarget  environment  the  primary  cause  of  perform¬ 
ance  degradation  of  the  existing  system  is  due  to  the  mismatch  between  the 
signal  processor  design  and  the  environment  in  which  the  signal  processor  must 
operate.  Therefore,  a  logical  approach  to  minimize  the  loss  of  performance  is 
to  examine  the  optimum  structure  of  the  signal  processor  under  a  multisensor, 
multitarget  environment.  Knowing  the  form  of  the  optimum  processor,  one  can 
then  explore  various  options  for  a  practical  system  realization. 


1.2  TECHNICAL  OBJECTIVES 


The  primary  objective  of  this  study  is  to  (1)  derive  the  optimum  time 
delay  estimator  under  a  multisensor,  multitarget  environment,  (2)  evaluate 
the  appropriate  performance  bound,  and  (3)  investigate  suboptimal  processor 
realizations . 


1.3  PREVIOUS  WORK 

Optimum  signal  processor  design  under  a  stationary  multisensor,  multi¬ 
target  environment  has  been  studied  by  a  number  of  researchers,  i.e.,  the  work 
by  Schweppe  (Reference  5)  on  sensor  array  data  processing  for  multiple  signal 
source,  Schultheiss  (Reference  6)  on  passive  sonar  detection  in  the  presence 
of  interference,  and  Anderson  and  Rudnick  (Reference  7)  on  rejection  of 
coherent  signal  arrival.  In  addition,  there  are  the  works  by  Capon  (Refer¬ 
ence  8),  Steinberg  (Reference  9),  Cox  (Reference  10),  McGarty  (Reference  11), 
Rockmore  and  Bershad  (Reference  12),  and  Owsley  and  Swope  (Reference  13). 

With  the  possible  exception  of  the  work  by  Schweppe  and  Ows 1 ey /Swope ,  the 
primary  efforts  of  these  studies  were  on  determining  the  effect  of  inter¬ 
ference  on  the  existing  (single  target)  processor.  Using  a  Least  Mean  Square 
approach,  Schweppe  derived  a  decoupled  beamformer.  However,  his  formulation 
ignored  the  effect  of  dependence  of  the  covariance  matrix  on  the  unknown 
parameter.  On  the  other  hand,  Owsley  and  Swope  derived  a  single  frequency 
multichannel  focused  beamformer  using  a  Weighted  Least  Square  approach. 
However,  Owsley  and  Swope's  formulation  did  not  include  the  bias  correction 
term. 


1.4  TECHNICAL  APPROACH  AND  ORGANIZATION 

This  study  provides  a  fundamental  examination  of  the  optimum  signal  pro¬ 
cessor  design  for  time  delay  estimation  under  the  assumption  of  a  multisensor, 
multitarget  environment.  Using  an  MLE  procedure,  an  optimum  multi  sensor, 
multitarget  time  delay  estimator  is  derived.  The  resulting  signal  processor 
is  reduced  to  its  simplest  form  for  system  realization.  In  addition,  this 
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study  derives  the  appropriate  performance  bound  for  the  resulting  estimator. 
Comparisons  between  optimum  and  suboptimum  realizations  are  also  discussed. 

The  organization  of  this  report  is  as  follows:  Chapter  2  discusses  and 
formulates  the  multisensor,  multitarget  time  delay  estimation  problem.  Chap¬ 
ter  3  derives  the  optimum  time  delay  estimator  and  an  appropriate  performance 
bound.  In  addition,  the  extension  of  time  delay  estimation  to  localization 
parameter  estimation  and  power  spectral  estimation  is  considered.  Chapter  4 
provides  some  discussion  on  suboptimum  processor  realization.  Chapter  5 
discusses  alternate  approaches  for  improved  multitarget  parameter  resolution. 
Chapter  6  presents  the  sunmary  and  the  conclusions  of  the  study. 


3/4 
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CHAPTER  2 

PROBLEM  FORMULATION 


2.1  INTRODUCTION 

In  this  chapter  we  discuss  the  mathematical  description  of  the  sensor 
observations  and  formulate  the  general  multisensor,  multitarget  time  delay 
estimation  problem.  The  optimum  signal  processor  is  sought  via  an  MLE  pro¬ 
cedure.  In  general,  we  assume  the  target/sensor  environment  consists  of 
M  sensors  and  J  targets  (or  J  arrival  paths).  Target  waveforms  are  assumed 
mutually  uncorrelated.  As  an  introduction  to  the  general  multitarget,  multi¬ 
sensor  problem  formulation,  we  first  consider  the  basic  description  of 
observables  for  a  single  target  case. 


2.2  DESCRIPTION  OF  OBSERVABLES 

Let  M  sensor  arrays  be  distributed  arbitrarily  in  space.  These  M  sensor 
arrays  then  produce  M  continuous  waveforms  which  represent  the  M  spatial 
samples  of  a  random  field  generated  by  the  target  signal  in  additive  ambient 
noise.  Let  the  M  waveforms  be  observed  for  a  duration  of  T  seconds  and  let 
the  received  waveforms  be  written  as: 


y(t,  2,.)  =  a.  s.(t)  +  n . ( t )  ;  i  =  1,  2 . M 

1  11  1  t  e  [0,  T] 


(2.2-1) 


where  a^ ,  s.(t),  and  n^t)  are  the  signal  attenuation  factor,  signal  waveform, 
and  the  noise  waveform,  respectively,  for  the  ith  sensor.  We  note  that  the 
observed  waveform  y(t,  a.)  is  a  random  function  of  space  and  time,  and  a.  is 
the  location  vector  of  the  ith  sensor  with  respect  to  a  known  reference  point. 
(See  Figure  2-1.) 
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Figure  2-1.  Description  of  Passive  Sensor  Observation 

To  simplify  the  analysis,  the  following  assumptions  have  been  made: 

1.  Pure  Time  Delay  Channel  -  The  signal  waveforms  received  at  each 

sensor  are  identical  except  for  a  pure  time  delay.  Thus,  the  signal 
waveform  (for  the  single  target  case)  of  the  ith  sensor  can  be 
written  as  s^(t)  =  s(t  +  D. ),  where  Di  is  the  propagation  delay  from 
the  signal  source  to  the  ith  sensor.  Therefore,  Equation  (2.2-1) 
can  be  rewritten  as: 


y(t,  H.)  =  a,  s ( t  +  D,)  +  n i ( t )  ;  i  -  1,  2,  ...,  M  (2.2-2) 
1  1  1  t  e  [0,  T]  . 


From  Figure  2-2,  we  see  that  the  propagation  time  delay  is  giv^n  by 
the  relation 


r  - 


c 


(2.2-3) 
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where 

r.  =  position  vector  of  the  signal  source 
£.j  =  location  vector  of  the  ith  sensor 
c  =  propagation  speed  of  the  medium. 

2.  Noise  at  each  sensor  is  assumed  additive  with  known  spatial  and 
time  correlation  function. 

3.  The  observation  time  T  is  much  longer  than  the  travel  time  of  the 
wavefront  across  the  sensor  array. 

4.  Both  the  signal  process  and  the  noise  process  are  mutually  uncor¬ 
related,  zero  mean,  Gaussian,  stationary  in  time,  homogenous  in 
space  (Yaglom  (Reference  14)),  and  have  a  known  band-limited  power 
spectrum. 

We  note  that  the  set  of  observations  y(t,  _?,.);  i  =1,  2,  ....  M; 
t  e  [0,  T]  contains  a  complete  description  of  the  observables.  However,  it  has 
infinite  dimensions  and  is  analytically  intractable.  Therefore,  we  wish  to 
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Figure  2-2.  Target -to-Sensor  Geometry 
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represent  the  time-limited  observation  of  the  waveform  by  a  finite  set  of  dis¬ 
crete  random  variables.  Because  of  the  time  delays  in  the  signal  waveform,  a 
convenient  approach  is  to  represent  the  observation  in  terms  of  a  Fourier 
series  expansion  as  discussed  by  MacDonald  and  Schultheiss  (Reference  15), 
Knapp  and  Carter  (Reference  1),  and  Rockmore  (Reference  16).  The  resulting 
coefficients  are  almost  uncorrelated  for  a  large  bandwidth  time  product  as 
demonstrated  by  Hodgkiss  and  Nolte  (Reference  17).  Therefore,  Equation 
(2.2-2)  can  be  represented  by  a  Fourier  coefficient  vector  as: 


2k 


k  =  ±1,  ±2, 


(2.2-4) 


where 


2k  =  (<*lk  ^k 


/  H°1  „  aja)kD2  „  eJVM Y 

^k  "  \  !  2e  aMe  / 


Ac  =  Wk  ^k  Alk^ 


(2.2-5a) 


(2.2-5b) 


(2.2-5c) 


are  the  vectors  of  observation,  signal  steering,  and  additive  noise,  respec¬ 
tively,  for  all  sensor  outputs  at  frequency  =  2-rrk /T.  Furthermore,  the  var¬ 
iables  a-^,  and  are  defined  as  follows: 


aik  “  T  /  dt 

Jo 


(2.2-6a) 
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-joi,1 

s(t)  e  dt 


1 

^ik  '  T 


ni(t) 


(2.2-6b) 


(2.2-6c) 


Note  that  for  real  signal  and  noise  waveforms,  the  Fourier  coefficients  are 
conjugate  symmetric;  i.e.,  _k  =  ct*k.  Thus,  for  a  band-limited  process  one 
only  needs  to  consider  frequency  components  k  =  1,  2,  ....  B,  where  B  is  the 
highest  cut-off  frequency  of  either  the  signal  or  the  noise  spectrum.  For  a 
low  pass  power  spectrum,  B  can  be  identified  as  the  one-sided  bandwidth. 

Since  both  the  signal  and  noise  are  zero  mean  Gaussian  processes,  it  is 
easily  established  that  is  a  zero  mean  Gaussian  random  vector  with  covari¬ 
ance  matrix: 


Rk  =  E(°k  s$) 


-  Sk  h  X  + 


Nk  ^k 


(2.2-7a) 


where  a£  is  the  complex  conjugate  transpose  of  a^;  Sk  and  Nk  are  the  dis¬ 
crete  signal  and  noise  power  spectral  density,  respectively,  at  frequency 
u^;  and  Qk  is  the  normalized  spatial  covariance  matrix  of  the  noise  vector 
such  that 


tr(Qk)  =  M  (2.2-7b) 

where  tr(  )  defines  the  trace  of  a  matrix. 

Now  writing 
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B 

=  ,  MBn  lRk  I _1  exP{-«k  Rk  “k}  (2.2-10) 

k=l 

as  the  pdf  of  the  complete  observation  vector  a  conditioned  on  the  propagation 
time  delay  vector  D.  Thus  Equation  (2.2-10)  provides  a  complete  statistical 
description  of  the  observation  conditioned  on  the  unknown  time  delay  vector. 

2.3  STATEMENT  OF  THE  PROBLEM 

In  this  section-  we  formulate  the  general  multisensor,  multitarget  time 
delay  estimation  problem  in  terms  of  the  observation  vector  described  in  the 
previous  section.  We  consider  the  problem  of  optimum  time  delay  estimation 
from  M  sensor  arrays  in  the  presence  of  J  possible  targets.  The  final  objec¬ 
tive  of  sensor  array  time  delay  estimation  is  to  provide  optimum  estimates  of 
target  localization  parameters  (i.e.,  range  and  bearing).  We  assume  that 
interarray  separation  is  large  compared  to  the  size  of  each  sensor  array. 

Such  an  arrangement  is  normally  found  in  large  aperture  surveillance  systems. 
See,  for  example,  the  large  seismic  array  described  by  Capon  (Reference  18). 

A  recent  study  by  Carter  (Reference  2)  has  shown  that  for  passive  localization 
of  an  acoustic  source,  the  variance  of  the  bearing  error  is  inversely  propor¬ 
tional  to  the  square  of  the  base  length  while  the  variance  of  the  range  error 
is  proportional  to  the  fourth-power  of  the  ratio  of  the  true  range  to  the  base 
length.  Thus,  a  long  base  length  is  desirable  in  reducing  the  variance  of  the 
estimates. 

Let  the  M  sensor  array  outputs  from  J  acoustic  sources  be  written  as: 

J 

y(t,  ii)  =23  aij  Sjft  +  Dij)  +  ni(t)  ;  t  e  [°,  T]  (2.3-1) 

j=l 
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where  i  =1,  2,  M;  a^-  is  the  signal  attenuation  factor  for  the  jth  tar¬ 
get  to  the  ith  sensor;  D^-  is  the  propagation  time  delay  from  target  j  to 
sensor  i;  and  s j ( t )  and  n-(t)  are  the  band-limited  signal  and  noise  processes 
with  the  usual  assumptions  of  zero  mean,  Gaussian,  time  stationary,  and  spati¬ 
ally  homogeneous.  Fourier  expansion  of  Equation  (2.3-1)  yields  the  finite 
dimensional  frequency  domain  representation  as: 


=  wiA  +  A<  ; 


k  =  1,  2,  ...,  B 


(2.3 -2a) 


where 


2k  =  (alk  “2k  “Mk^ 


(2.3-2b) 


4  =  (Bkl  Sk2  8kj) 
3c  =  (nlk  ^k  •**  ’Vlk^ 


(2.3-2c) 


( 2 . 3 -2d ) 


i  =1,  2,  . . . ,  M 
j  =  2,  ...,  J 


(2.3-2e) 


is  a  M  x  J  complex  "propagation  delay"  matrix,  where  M  is  the  number  of 
sensors,  J  is  the  nuntoer  of  targets,  and  =  2irk/T  is  the  discrete  frequency. 
Similar  to  the  development  in  the  last  section,  the  observation  vector  is 
Gaussian  with  the  following  statistics: 


E(cfc)  =  0 


(2.3-3a) 


E(2k  2k)  = 


C;  if  k  t  l 
;  if  k  =  l 


( 2 . 3 -3b ) 


TR  6757 


where 


Rk  A  Wk  Sk  W£  +  Qk 


(2.3-3c) 


and 


Sk  =  E(£k  * 

For  uncorrelated  sources,  we  have 


(2.3-3d) 


S|<  =  d i ag { Sk ^ ,  Sk2,  ...  ,  SkJ} 


(2.3-3e) 


as  a  J  x  J  diagonal  matrix 
of  target  j.  Thus,  the 
defined  by 


with  element  Skj  as  the  signal  power  at  frequency 
pdf  of  conditioned  on  the  time  delay  matrix 


D  =  {0_j  j }  ;  i  -  1,  2,  ....  M  (2.3-3f ) 

j=l,  2,  ...,  J 
is 

p(°!k  I D)  "  \\\~l  exp{-af  R,;1  ^ }  .  (2.3-4) 

Finally,  writing 

a  =  (aj  £  ...  aJjT  ,  (2.3-5a) 


then  the  pdf  of  a,  the  complete  observation  vector,  conditioned  on  D,  the  mul¬ 
tisensor,  multitarget  time  delay  matrix,  can  be  written  as 


B 

p(ot|  D)  n  p(ak |D) 
k=l 

B 

=  ,  MBn  |Rk!_1  exp{-a£  R'1  .  (2.3-5b) 

k=l 
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Equation  (2.3-5b)  constitutes  the  basis  for  the  derivation  of  the  multisensor, 
multitarget  optimum  time  delay  processor.  An  MLE  is  obtained  by  maximizing 
the  conditional  pdf  with  respect  to  each  element  of  the  propagation  time  delay 
matrix.  Thus,  symbolically  one  can  write 

=  Max  Arg  p(a|D)  .  (2.3-6) 


14 


TR  6757 


CHAPTER  3 

OPTIMUM  MULTITARGET  PARAMETER  ESTIMATION 


3.1  INTRODUCTION 

In  principle  the  optimum  estimate  of  the  time  delay  matrix  D  is  obtained 
by  solving  Equation  (2.3-6)  numerically.  In  general,  a  direct  implementation 
of  Equation  (2.3-6)  results  in  a  very  complex  processor.  However,  without 
loss  of  performance,  Equation  (2.3-6)  can  be  reduced  to  its  simplest  form  by 
mathematical  manipulations.  The  resulting  processor  is  usually  realizable. 

In  this  chapter  we  investigate  the  fine  structure  of  the  optimum  multi¬ 
target  time  delay  signal  processor.  Since  time  delays  are  modeled  as  unknown 
constants,  we  seek  an  optimum  estimator  via  an  MLE  approach.  We  provide  a 
performance  bound  for  the  resulting  estimator.  We  study  in  detail  the  two¬ 
sensor,  two-target  case.  For  the  three-sensor  case  we  establish  the  relation¬ 
ship  between  optimum  localization  parameter  estimation  and  optimum  time  delay 
estimation.  Finally,  we  investgate  the  structure  of  the  optimum  power  spec¬ 
tral  estimator  when  the  target  signal  power  spectrum  is  not  known  a  priori. 


3.2  THE  LIKELIHOOD  EQUATION 

It  was  shown  in  Chapter  2,  Equation  (2.3-5b),  that  the  pdf  of  the  sensor 
observation  vector  &  conditioned  on  the  time  delay  matrix  D  is  given  by 

B 

p(a|D)  n  p(ak!D)  (3.2-la) 

k=l 


where 


p(ok|D)  =  ir-M|Rkr1  exp{-o£  R^1  o^}  . 


(3.2-lb) 
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The  ij  element  of  the  time  delay  matrix  D  is  given  by  (see 
Equation  ( 2 . 3-3f ) ) 


i  =  1,  2,  ....  M 
j  ~  1 »  2 ,  « . .  >  0 


(3.2-2) 


where  r.  and  denote  the  vector  locations  of  target  j  and  sensor  i,  respec¬ 
tively.  The  log-likelihood  function  of  Equation  (3.2-la)  is  defined  by 


B 

A(D)  =  X]  1og  p(-k 
k=l 

B 

=  -MB  log(ir)  -^Ak(D) 
k=l 


where 


VD>  =  log|Rk| 


(3.2-3a) 


(3.2-3b) 


and  Rk  is  the  covariance  matrix  of  the  zero  mean  complex  Gaussian  vector 
given  by  (see  Equation  (2.3-3c)) 

Rk  =  Wk  Sk  Wk  +  Nk  Qk  (3.2-3c) 

where  for  the  uncorrelated  source  case, 

Sk  =  diag{Skl,  Sk2,  ...,  SkJ>  (3.2-3d) 


is  a  J  x  J  diagonal  signal  power  density  matrix  whose  jth  diagonal  entry 
denotes  the  discrete  signal  power  density  of  target  j  at  frequency  tok.  Nk  is 


the  discrete  noise  power  density,  and  is  the  normalized  covariance  matrix. 
Finally,  the  M  x  J  "propagation  delay  matrix",  W^,  is  defined  by 


i  =  1,  2,  ....  M 

j  —  1 »  2,  ...,  J 


(3.2-3e) 


where  a^  denotes  the  signal  attenuation  factor  from  target  j  to  sensor  i. 
Now  writing 


“k-«U  *k2-w 


(3.2-4a) 


where 


(v 


nDu 


a2je 


J“k°2j 


aMje 


k  Mj 


(3.2-4b) 


is  the  attenuated  signal  steering  vector  of  target  j  at  frequency  u^.  Using 
the  fact  that  is  a  diagonal  matrix  for  uncorrelated  sources,  we  have  the 
obvious  identity: 


wksk«i;  =  Eskj  v»k*j 

j=l 


0 

l 

j=l 


Skj  Pkj 


(3.2-4c) 


where  .  is  defined  by 
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Note  that  we  have  chosen  to  let 


(3.2-4d) 


or  equivalently  (to  reference  all  time  delays  to  the  first  sensor),  one  can 
wr  i  te 


=  (alj 


...  aMje 


HVl.j 


(3.2-4e) 


where  A.  ,  is  given  by  the  relation 

*  J 


Aij  =  °i+l, j  "  Dlj  ;  1  -  1,  2,  ....  M-l  . 


(3.2-4f) 


Now  .  is  an  M  x  M  Hermitian  matrix  of  rank  one  and  is  a  function  of  the 

attenuation  vector  a.  and  the  time  delay  vector  A..  The  vectors  a_.  and  A-  are 

J  J  J  J 

given  by 


-j  =  (alj  a2j  aMj} 


(3.2-4g) 


-j  =  {Alj 


(3.2-4h) 


Note  that  given  M  sensors,  there  are  M-l  independent  time  delay  pairs  from  a  pos¬ 
sible  total  of  M(M-l)/2.  The  selection  of  this  set  is  not  unique.  However,  for 
time  delay  estimation,  it  is  reasonable  to  assume  a  set  with  minimum  total  delay. 
This  set  (for  a  line  array)  is  given  by  the  inter-sensor  time  delay  vector: 

-j  =  (xlj  x2j  *•’  xM-l,j)T  ;  j=1’  2’  J  (3.2-4i ) 

where 
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Note  that  the  A  • • ‘ s  can  be  expressed  in  terms  of  the  t- .'s.  For  example: 

*  J  '  J 


Aij  =  "  °lj 


~  Dn+l,j  "  Dnj 


=  uT  x. 

-i  -j 


(3.2-4J) 


where  is  a  column  vector  whose  first  i  entries  are  one  and  the  remainder 
are  zero. 

Furthermore,  by  writing  the  mn  element  of  the  matrix  Pkj  as  P™,  and 
from  Equation  (3.2-4e),  we  can  identify  the  relation 


Pmn  -  a  a  j'Vl.j)  . 

Kkj  amj  anj  e 


m  =  1,  2,  M  (3.2-5a) 

n  —  X,  2,  M  . 


Now  utilizing  the  relation  in  Equation  (3.2-4j),  we  obtain 


Pmn  =  a  a  -J 

Kkj  amj  anj  e 


(3.2-5b) 


Thus  we  have  shown  that  the  observation  covariance  matrix  can  be  expressed 
as  a  function  of  inter-sensor  time  delays  instead  of  the  actual  propagation 
delays.  Therefore,  the  direct  observable  quantities  are  time  delays  instead 
of  the  propagation  delays. 


Using  Equation  (3.2-4c)  in  (3.2-3c),  one  obtains  an  alternate  expression 
for  the  observation  covariance  matrix: 
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\ =  E  skj  pkj  *  \  %  ■ 
j=i 


(3.2-6) 


Equation  (3.2-6)  can  also  be  summed  in  the  following  way: 


Rk  =  Skl  Pki  +  Nkl  Qkl 


=  Sk2  Pk2  +  Nk2  Qk2 


~  SkJ  PkJ  +  NkJ  QkJ 


(3.2-7a) 


where 


*kj  =  ^2  Ski  +  \  ;  3  *  1,  2,  ....  J 

i-1 

1« 


(3.2-7b) 


2  Ski  Pki  Nk  % 

i=l 

i^j 


/Nkj 


(3.2-7c) 


Note  that  Qy  can  be  regarded  as  an  equivalent  noise  process  and  is  indepen¬ 
dent  of  the  parameter  t. •  For  simplicity  we  shall  assume  the  propagation 

J 

attenuation  coefficients  are  known  and,  for  convenience,  assume  they  are  equal 
to  unity. 

Finally,  writing  the  incremental  (intersensor)  time  delay  vector  for  all 
targets  as 


e  s  (tT  tT  .TjT 

0  VTj  Tg  ...  TjJ  , 


(3.2-8) 
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a  p-parameter  column  vector  where  p  =  J(M-l),  the  vector  log-likelihood 
function  (Equations  (3.2-3a)  and  (3.2-3b))  becomes 

B 

A(e)  =  -MB  log(ir)  -  A^e)  (3.2-9) 

k=l 

where 

\(Q)  =  log  |Rk  |  +  c£  Rk_1  .  (3.2-10) 

The  MLE  of  e  is  one  that  maximizes  the  likelihood  function  A(e).  A  neces¬ 
sary  condition  for  the  location  of  the  maximum  is  given  by  the  vector  likeli¬ 
hood  equation 


B 

VA(e)  =  VA^e)  =  0  (3.2-11) 

■  k=l 


where  7  is  the  gradient  operator  given  by  the  column  vector 

7*  1i^)T  •  <3'2-12> 

3.3  MULTIPLE  PARAMETER  ESTIMATION 

It  was  shown  in  the  previous  section  that  the  optimal  estimate  of  the 
time  delay  vector  is  hinged  on  solving  the  vector  likelihood  equation.  For 
the  MLE,  we  have 


B 

£(»)  =  VA(©)  =  -  ^  VAk(e)  =  0  (3.3-1) 

k=l 


where  f(e)  is  a  vector  function.  The  likelihood  function  Ak(e)  is  defined  as 
in  Equation  (3.2-10). 
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In  this  section  we  examine  the  solution  of  the  likelihood  equation  in 
general  and  in  Section  3.5  we  explore  in  detail  the  structure  of  the  optimum 
multitarget,  multisensor  processor. 

Let  e j  denote  the  jth  element  of  the  time  delay  vector  e,  then 
Equation  (3.3-1)  can  be  written  as 


_3 

ae 


u 

:  Me)  =  af-  ^(e)  =  0  ; 


j  ”1,  2,  ...,  p 


(3.3-2) 


k=l 


From  Equation  (3.2-10)  we  obtain 


3Rk 

+  air—  °5k 

J 


3r:a 

J 


(3.3-3) 


where  tr(  )  denotes  the  trace  of  a  matrix.  Note  that  the  evaluation  of  the 
first  term  is  a  straightforward  application  of  the  chain  rule.  The  exact  pro¬ 
cedure  can  be  found  in  Rockmore  (Reference  16).  The  derivative  of 
Equation  (3.3-3)  w.r.t.  ei  is 


3% 

aeTaeT  v- 

'  J 


(e)  =  tr 


( ^k1 3Rk  +  .-i  »\  \ 

\3e.  aej.  Kk  3ei3e,  / 


a2^1 

+  ^k  3eTaeT  2k 


(3.3-4) 


Using  the  linear  property  of  the  trace  operator  and  the  relation 


3R, 


-1 


ae. 


iR.l 

3e.  Rk 


» 
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Equation  (3.3-4)  becomes 


32Ak(e) 

38.38. 

J 


=  -tr 


/  -1  3Rk  -1  9Rk  \ 

(R* 


+  tr  R 


-1 


*\ 

3e.3e . 

*  J 


a2^1 

+  ^k  3e  .3e .  ^k 

'  «J 


(3.3-5) 

Equation  (3.3-5)  is  important  in  evaluating  the  performance  bound  of  an  esti¬ 
mator.  Substituting  Equation  (3.3-3)  in  (Equation  3.3-1),  one  obtains  the 
set  of  necessary  conditions  for  the  ML  estimate: 


fi(e)  = 


k=l  '  1  '  1 


=  0  ;  i  =  1,  2,  ...,  p. 

(3.3-6) 


Equation  (3.3-6)  is  usually  non-linear.  The  structure  of  the  optimum 
processor  can  be  found  by  reducing  the  required  mathematical  operations  to 
simplest  form.  In  principle.  Equation  (3.3-6)  can  be  solved  by  searching  the 
p-parameter  space  for  a  simultaneous  null.  A  more  efficient  algorithm, 
however,  is  implementing  a  closed-loop  null  tracker.  Further  discussion  on 
this  important  subject  is  beyond  the  scope  of  this  study. 


We  remark  that  Equation  (3.3-6)  is  the  necessary  condition  for  the 
existence  of  a  maximum.  For  sufficiency  it  requires  not  only  the  condition 
given  in  Equation  (3.3-6)  but  also  the  following  (Bryson  and  Ho, 

Reference  19): 


L_  f(s)  =  3  A(e)  g 
3e  US)  9q2  ‘ 


(3.3-7) 


2  2 

By  this  we  mean  the  square  matrix  3  A(e)/3e  must  be  negative  definite. 
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3.4  ESTIMATOR  PERFORMANCE  EVALUATION 

In  this  section  we  derive  the  multi -parameter  Cramer-Rao  Lower  Bound 
(CRLB)  and  show  that  the  resulting  estimate  obtained  from  solving 
Equation  (3.3-6)  satisfies  the  bound  for  a  large  observation  time.  There¬ 
fore,  the  resulting  estimate  is  efficient. 

The  CRLB  for  an  unbiased  estimate  of  the  ith  parameter  is  given  by  Van 
Trees  (Reference  20): 


(3.4-1) 


where  (  )..  denotes  the  ith  diagonal  element  of  a  matrix  and  J  is  the  Fisher 
Information  Matrix  whose  ij  element  is  defined  by 


=  E  (im 


i  =  1,  2,  ....  P 
j  ~  2,  ...,  P 


(3.4-2) 


where  A(e)  is  the  log-likelihood  function.  Using  Equations  (3.3-2)  and 
(3.3-5)  in  (3.4-2),  we  immediately  obtain 


A  /  j  3Rk  j  3Rj<  \  /  _i  9  \  ^ 

iJ  =  -2J  tr(Virv4)-tr  Virs7  +  scir 

'  1  J  /  '  1  J  1  J 


i  =  -eo 


(3.4-3) 


where  ©  denotes  the  true  parameter  value.  Equation  (3.4-3)  can  be  simplified 
as  follows:  Taking  the  derivative  of  the  identity  Rk  R)<=I  first  w.r.t.  ©j 
and  then  w.r.t.  ©-,  one  obtains  the  relation 


,  32R.  32R^ 

q-l _ + _ !L_  r 

\  3©  •  3©  ■  3©  •  3©  •  k 

I  J  J 


3R^  3Rk  SR^1  3Rk 

le .  3e~  3© .  3© . 

J  l  *  J 


(3.4-4) 
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Substituting  Equation  (3.4-4)  in  Equation  (3.4-3),  the  latter  can  be  simpli¬ 
fied  to 


(3.4-5a) 


(3.4-5b) 


Note  that  because  tr(AB)  =  tr(BA),  we  have  J. .  =  J...  Equation  (3.4-5b) 

'  J  J  ■ 

agrees  with  the  expression  obtained  by  Bangs  (Reference  21).  However,  the 
derivation  presented  here  is  somewhat  simpler  and  more  direct. 

A 

Next  we  proceed  to  show  that  the  MLE,  e^L,  of  e  obtained  by  simultan¬ 
eously  solving  the  set  of  equations  in  Equation  (3.3-6)  is  unbiased  and 
achieves  the  CRLB. 

Let  f(e)  *  (fx(e)  f2(e)  ...  fp(e))T 

and  write  the  Taylor  series  expansion  of  f(e)  about  the  true  parameter  value 

as  follows: 

— o 


!<!>  *  i<»0>  *  ^ 


"  i)0  + 


9  =  0_ 
—  -0 


(3.4-6) 


A 

Since  f(e^)  =  0  by  definition,  we  have  after  neglecting  higher  order  terms 
(small  random  error  assumption): 


3KV  - 

3e  ^®ML  "  2o)  * 


(3.4-7) 
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Now  assume  the  law  of  large  number  applies  so  that  (for  a  sufficiently 
long  observation  time)  one  can  replace  the  derivative  by  its  expected  value. 
Taking  the  expected  value  on  both  sides  of  Equation  (3.4-7),  one  obtains 


“  J  Ef  3e~° 


3f(ejl  )  -1 


E[f (Sq)] 


(3.4-8) 


But  from  Equation  (3.3-6) 


=  -E  trklw:)  +  tr(S~ 

k=l  L  1  v  1 


[  /  _i  9\  9Ri<1  \l 

-L,  tr  y v  air +  r-  Rk )  fl_e 
k=i  L  \  i  i  /  J  !-?o 


*  0  . 


(3.4-9) 


Therefore,  ^(e^^)  =  e^,  and  is  an  unbiased  estimator.  Post  multipling 
Equation  (3.4-7)  by  its  conjugate  transpose  and  taking  the  expectation  yields 


3f(eJ 


— O'—  '-O' 


r-  *  t  3f  (s^1 

E  r^ML  "  V  "  SqM  E  9e~ 


(3.4-10) 


Now  recall  that  £(e)  =  VA(e),  so  one  can  write 


32A(e) 


i  =  ^o 
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32  A(  © ) 

32A(e) 

32A(@) 

3ex2 

3©^3©2 

30  ,  30 „ 
1  P 

32A( © ) 

32A(0) 

32A(0) 

38239^ 

3©  2 

"  3e239p 

32A(e) 

Vel 


32A(©) 

3e„3e9 
P  2 


Furthermore,  we  have  the  relation 

E[f(9o)  fM^)]  «  E[7A(e)  VTA(9)] 


£  =  ®0 


=  J  . 


32A(o) 


(3.4-11) 


9 


(3.4-12) 


Therefore,  substituting  Equations  (3.4-11)  and  (3.4-12)  in  (3.4-10)  yields 
E  "  9o^  =  J  J 


(3.4-13a) 


VAR(9i)ML  - 


( 3 .4 -13b) 


Thus  we  have  shown  that  the  multi -parameter  MLE  is  an  efficient,  unbiased 
estimator  for  large  observation  time. 
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yi 
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3.5  MULTISENSOR,  MULTITARGET  PARAMETER  ESTIMATION 

In  this  section  we  investigate  the  detailed  structure  of  the  optimum 
signal  processor  for  multiple  parameter  estimation.  For  the  purpose  of  target 
localization,  the  parameter  set  of  interest  is  time  dr^ys,  range,  and 
bearing.  There  are  two  possible  approaches  in  estimating  the  localization 
parameters.  The  first  is  via  a  geometric  transformation  from  measured  time 
delays;  the  second  is  via  a  direct  range  and  bearing  signal  processor.  From 
the  results  of  this  study,  we  will  clarify  the  relative  merits  between  these 
two  approaches.  In  Section  3.5.1,  we  derive  the  optimum  multisensor,  multi¬ 
target  time  delay  processor.  In  Section  3.5.2  we  discuss  the  various  method¬ 
ologies  of  obtaining  range  and  bearing.  Finally,  in  Section  3.5.3  we  briefly 
discuss  the  problem  of  optimum  time  delay  estimation  with  unknown  target  power 
spectra. 


3.5.1  Time  Delay  Estimation 

It  was  shown  in  Section  3.3  that  the  MLE  of  the  time  delay  parameter 
vector  e  required  the  simultaneous  solution  of  the  vector  likelihood  equation: 


g 

^  =0;  j  -  1,  2,  ....  J(M  -  1)  (3.5.1-la) 

k=l 


where 


9R, 


■1 


3Ak(-}  _  (p-l  9Rk  \  *  Zt- 

3e .  '  tr  I  Pk  3e .  )  ^k  3e  .  • 

J  '  J  /  J 


(3.5.1-lb) 


We  will  explore  the  detailed  structure  of  the  optimum  signal  processor  by 
simplifying  Equation  (3.5.1-la).  Recall  from  Equation  (3.2-6)  that 


Rk  =  Skj  Pkj  +  Nkj  Qkj 


(3.5.1-2a) 
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T5 


N 


b 


where 


qkj 


Ski  Pki  +  Nk  Qk 


/Nkj 


(3.5.1 -2b) 


and 


"kj  ■  E  ski +  \ 

i=l 


(3.5.1 -2c ) 


For  notational  simplicity  we  shall  assume  in  Equation  (3.5.1-la)  that  there 
are  J  targets  but  with  two  sensors.  The  e.'s  in  this  case  correspond  to 

J 

the  time  delay  for  each  target.  In  the  case  of  more  than  two  sensors,  9j 
must  be  replaced  by  each  element  of  the  inter-sensor  time  delay  vector 
(t j j ,  T2j,  . ...  tM-1  j)  of  target  3.  where  t^-  denotes  the  time  delay  between 
sensors  i  +  1  and  i  of  target  j. 

Applying  the  well-known  matrix  inversion  Lemma  (HTR‘^H  +  = 

M  -  MHT(HMHT  +  R)_1HM  to  R^1,  we  obtain 

"k1  ■  <Skj  pkj  *  \i  Qkj  >_1 


\j  ltGkjVNkj 


qkj  pkj  Gkj 


(3. 5. 1-3) 


«-l 


where  G^.  =  Q^.  is  defined  as  the  effective  array  gain  for  target  j. 

The  derivative  of  w.r.t.  9j  can  be  written  as 


3R 

F. 

3e. 


-1 


/!!ni 

\  ‘ 


akj  Wf  Pkj)  Qkj 


(3.5.1~4a) 
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where  we  have  defined 


K-i 


kj 


W?1 


1  *  Gkj  sk0/Nkj 


(3.5.1 -4b  J 


akj  "  V,>k/  • 


(3. 5.1 -4c) 


Using  Equations  (3. 5. 1-3)  and  (3.5.1-4a),  the  first  term  in  Equation 
(3.5.1-lb)  can  be  written  as  follows: 


--1  3?h 

k  '  "kj  '  '  kJ  3eJ 


-  skj|hkjl2  tr  (Gk1  pkj  Gk1  If 


(3.5.1-5a) 


But  this  can  be  further  simplified  using  the  following  relations: 


tr /g-1  !!k\,  tr  I" q-1  ( v*  +v  *3l\ 

38j  )  tr  qkj  \  3e.  hs  +  hi  wr) 


4*j  qkj  +  air1  qkj  hi 

J  J 


aGkj 

30j 


(3.5.1-5b) 


and 


tr  («}  pkj  «k1 £f )  '  *'  [«k1  4*0  %'j  (If  4*0  *  4o  If ) 


3G.  • 
=  G,  .  — ^ 


ki  a«.  * 


(3.5. 1 -5c ) 
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Thus  Equation  (3.5.1-5a)  reduces  to 


3R 


3Gi, 


tr  r; 


1  \  _  kj  oukj 


39j 


Nkj  38j 


’kj  l"kj 


|2  5  fki 

i1  bkj  3e . 


(3. 5. 1-6) 


Using  Equations  (3.5.1-lb),  (3.5.1-4a)  and  (3.5.1-5c)  in  (3.5.1-la),  the 
likelihood  equation  becomes 


B 

3A(e)  _ 

3e .  "  Z-rf 

J  k«l 

=  0. 


kj 


$ 


3Gk; 


akj  3ejJ  Pkj)  Gkj 


j  1  •••  J 


-  3G|< 

*kj  9^ 


(3. 5. 1-7) 


Equation  (3. 5. 1-7)  reduces  to  that  obtained  by  Bangs  (References  21,  23)  for 
the  single  target  case.  It  should  be  pointed  out  that  our  development  up  to 
here  in  many  ways  parallels  Bang's  work.  However,  there  are  also  major 
differences.  We  are  interested  in  a  multitarget  environment  while  Bangs' 
work  dealt  exclusively  with  single  target.  We  are  interested  in  joint  time 
delay  vector  estimation  while  Bangs'  work  is  primarily  concerned  with  range 
and  bearing  estimation.  Thus,  our  work  in  this  section  can  be  considered  as 
an  extension  of  Bangs'  original  work  to  include  the  multitarget,  multisensor 
environment. 


For  time  delay  estimation,  the  likelihood  equation  (Equation  (3. 5. 1-7)) 
can  be  further  simplified  as  follows. 

Recall  from  Equations  (3.2-4e)  and  (3.2-4j)  that  the  steering  vector  for 
target  j  is 

v  .  .  (i  eJ“k“2IJ  ...  _  (3.5.1-8*) 

k  3 

where  the  attenuation  coefficients  a.,  were  assumed  known  and  for  convenience 

*  j 

they  were  assumed  to  have  unity. 
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Therefore,  the  mn  element  of  Pkj.  is 


m  ,  on(Vi->!n.l)TIj 


P  11  =  e 
*kj  e 


(3.5.1-8b) 


•HlVrVi’V 


•  -l  mn  ninri 

Ja,k$i  pkj 


(3.5.1-8c) 


where  $mn  =  (U,  ,  -  IT  t.  and  $rfn  is  defined  as  the  mn  element  of  the 
— m-i  — n-i  — j  i 

matrix  given  by 


„„  ..mn 
.mn  _  a$ 

1 


T  ^Ti 

(iU_i  -  u._i)  ; 


v*n-l  ^1-1'  5t77 

'  J 


i  =1,  2,  ....  M-l 
m  1,  2,  •  •  • ,  M 
n  =  1,  2,  ...»  M 


II  ;  if  n  <  i  <  m-l 
-1  ;  if  m  <  i  <  n-1 
0  ;  otherwise  . 


(3.5. l-8d) 


Note  that  for  a  system  of  M  sensors,  4^  is  an  M  x  M  matrix, 
let  M  =  3,  then 


For  example. 


0-1-1  0  0-1 
♦j  =  1  0  0  ;  $2=00-l. 

10  0  110 
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Define  1^  as  an  M  x  M  square  matrix  of  ones,  and  .  as  a  diagonal  "steering 
matrix",  whose  diagonal  elements  correspond  to  the  elements  of  the  steering 
vector 

Vkj  =  diag  j  1,  e  J  ,  e  J . e  1  J  j  .  (3. 5. 1-9) 


Then  the  following  relations  can  be  easily  established: 


P  =  V  IV* 
kj  vkj  vkj 


(3.5.1-10a) 


37^  ’  H  vkj  »t  vkj 


(3.5. 1 -10b ) 


Furthermore,  using  Equation  (3.5.1-10b),  Equation  (3.5. l-5b )  can  be  written  as 


3G.  .  ~  i 

w^; =  n  tr(Qkj  vkj  $i  vkj}  • 

*  j 


(3.5.1-10c) 


Finally,  using  Equations  (3.5.1-10a),  (3.5.1-10b),  (3.5.1-10c)  and 
replacing  e.  by  t-  .,  the  likelihood  equation  (Equation  (3. 5. 1-7))  can  be 

J  '  J 

simplified  to 


3A(t) 


fikj|Z  2k*  «k'j  VkjK  -  bT  !N)  %  °k‘j  2k  -  b) 


r1  a.  -  bkJ 


(3.5.1-lla) 


where 


blJ  *  *«  tr«kj  vkj  »i  vkj> 


(3.5.1-llb) 
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is  the  bias  correction  term  and  for  i  =  1,  2,  . ...  M-l  and  j  =  1,  2,  J 
where  M  is  the  number  of  sensors  and  J  is  the  number  of  targets.  Note  that  the 
J(M  -  1)  equations  are  coupled  and  must  be  solved  jointly  for  the  stationary 
point.  For  long  observation  time  such  that  the  frequency  samples  are  dense  over 
the  frequency  bands  of  the  signal  and  the  noise,  the  summation  in 
Equation  (3.5.1-lla)  can  be  replaced  by  integration  as  follows: 


3A(t) 
'  3?— 


CO  / 

^  vj(“)  $i  ■  bi(a)) 1 


V|(o))  Qj1  (to)  ot(a))/T  -  b'j(o))  j  dco  ; 

1  »  1,  2 . M-l 

j  -  2,  ...,  J. 


=  0  (3.5.1-12a) 

Note  that  in  obtaining  Equation  (3.5.1-12a)  the  relations  Tot.  =  a(u))  and 
| hk  |  =  T|h(u))p  have  been  used. 

It  is  recalled  that  for  the  jth  target 


S.(u>)/N^(uj) 


|  h .  (a>)  \  .  ~  - i 

J  1  +  G  .(a))  S, 


j(u)  S j  (a>)  /N j  (a») 


a  j  (ti>)  =  N  j  (to )  |hj(oj)  |‘ 


(3.5. l-12b ) 


(3.5.1 - 12c ) 


aj(o))  tr[Q^(to)  Vj(u))  ^(u)) 


VJ(o))] 


(3.5. 1 -12d ) 
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qj'h 


Nj(o)) 


U 

2  V“) 


N(u>)  Q ( oj) 


i=l 

L  i£j 


-1 


V“) 


J 

=  '^2  s-jM  +  N(“) 
i=l 
i^j 


Pj(«)  «  Vj(to)  1M  VJ(u) 
and 

6j(»)  -  VJ(«)  Qj1  ( u>)  Vj(id)  . 


(3.5.1 -12e ) 


(3.5.1 -12f ) 


(3.5.1 -12g ) 


(3.5.1 -12h ) 


The  optimum  multitarget,  multisensor  time  delay  signal  processor  is 
shown  in  Figure  3-1.  There  is  a  total  of  J(M  -  1)  processing  channels  for  the 
case  of  M  sensors  and  J  targets.  For  simplicity  we  show  a  single  processing 
channel.  Note  that  the  processing  channels  are  tightly  coupled.  The  signal 
conditioning  filters  depend  on  the  time  delay  parameters  from  other  process¬ 
ing  channels  as  well.  This  is  an  order  of  magnitude  more  complex  compared  to 
a  single  target  case.  A  number  of  suboptimal  realizations  can  be  found  as 
discussed  in  Section  4.  Finally,  we  remark  that  for  convenience,  we  show  the 
opti.ium  processor  (Figure  3-1)  in  the  continuous  frequency  domain.  For  prac¬ 
tical  considerations,  the  discrete  counterpart.  Equation  (3.5.1-lla),  is 
normally  used  since  the  correlation  process  can  be  mechanized  easily  via  the 
Fast  Fourier  Transform  (FFT)  algorithm. 

Using  Equations  (3.5.1 -12b-f ),  we  shall  study  the  specific  structure  of 
the  optimum  time  delay  processor  for  a  number  of  simple  but  important  cases. 


3. 5. 1.1  Case  1:  One  Target  and  Two  Sensors  (J  =  1,  M  =  2).  For  con¬ 
venience  we  assume  Q(co)  =  I;  i.e.,  the  noise  processes  are  equal  in  power  and 
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uncorrelated  between  sensors.  Noise  processes  being  uncorrelated  between 
sensors  is  a  reasonable  assumption  since  in  practice  sensors  are  separated  at 
least  at  a  half  wavelength  spacing. 


The  steering  vector  is  v 
fied  easily: 


Q”1  (w)  =  Q"1(oj)  =  I 


v_  =  (1  eJWT)^".  The  following  relations  can  be  veri 


(3.5.1-13a) 


G(cjo)  =  v*  Q"1  v  =  2 


(3.5.1 -13b) 


b(oj)  =  tr(Q_1  V  V*)  =  0 


(3.5.1 -13c ) 


|h(a)) 


S(oj)/N  a)) 

+  2  S(o)  /N(u) 


(3.5.1 -13d ) 


0-1  10 

;  V  * 

1  0  0  eJWT 


Thus  from  Equation  (3.5.1-12a),  the  likelihood  equation  is 


00 

Z(t)  =  j“|h|^  a*  V  V*  a  da)  =  0  . 


(3.5.1-13e) 


(3.5.1 -14a) 


This  simple  time  delay  processor  is  diagrammed  in  Figure  3-2.  This 
processor  is  identical  to  the  one  studied  by  Carter  (Reference  22)  and  can  be 
shown  as  follows.  From  Equation  (3.5.1 -14a) ,  we  have 

1  C  -  O  [1  0  1  [o  -ll  [l  0 

Z(t)  -  ir—  I  jculhl  ot*  .  .  a  daj 

^  l  "  Lo  eJ“TJ  Ll  0  J  Lo  e-0“TJ- 
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T 

=  H 


00 


|~|2  “l^  ofM 


T 

=  TH 


jui|h  |2  G^M  eJ0JT  dio 


-Jhf  !^l2  g12  ("I  eJ“XH 

-oo 


=  T 


3_ 

3t 


12 


(x) 


(3.5. 1 -14b) 


where  is  the  estimated  cross  power  spectrum  between  channels  a^(uj)  and 
<v,( u>),  and  Ri2(t)  is  the  GCC  function  studied  by  Knapp  and  Carter  (Refer¬ 
ence  1).  A  block  diagram  of  this  processor  is  shown  in  Figure  3-2 ( b ) .  Note 
that  the  null  of  Z(t)  corresponds  to  the  peak  of  R^2(t).  Furthermore,  the 
discrete  form  of  Equation  (3.5.1-14b)  is 

8 

Z(t)  «  2  Hlhk!2  «Lk  °2k  ^  (3.5.1  -14c ) 

k=-B 

B 

•IfE  lhklSk“?ke  V'  (3.5.1-14C) 

k=-B 

From  Equation  (3.5.1 -14b) ,  we  note  that  the  GCC  function  is  directly 
proportional  to  the  likelihood  function.  The  optimum  estimate  is  determined 
by  locating  the  peak  of  the  GCC  function  or,  equivalently,  the  null  of  its 
derivative.  We  have  shown  in  Section  3.4  that  the  MLE  is  efficient  for  a  long 
observation  time.  We  have  also  obtained  a  general  closed  form  expression  for 
the  CRLB.  The  CRLB  for  this  case  can  be  determined  easily.  From  Equa¬ 
tions  (3.4-1)  and  (3.4-5)  we  have 
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I 


■  2“kl\l2  sk 

and  the  CRLB  is 


VAR(t)  > 


(3.5.1-16a) 


(3.5.1 -16b) 


which  is  identical  to  the  expression  derived  by  Carter  (Reference  2).  Note 
that  T  is  the  observation  time  and 
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c(“>-(3ot¥uW)2  (3.6.1  -16c) 

is  known  as  the  magnitude  square  coherence  (MSC)  function  for  the  case  of  two 
equal-noise  power  channels. 


3. 5. 1.2  Case  2:  Two  Targets  and  Two  Sensors  (J  =  2,  M  =  2).  Again 
assuming  =  I  for  simplicity.  The  parameter  vector  consists  of  two 
elements;  i.e.,  e  =  (t-p  t2)^,  the  time  delays  to  target  number  one  and  target 
number  two.  The  optimum  estimates  can  be  obtained  by  solving  simultaneously 
the  two  likelihood  equations: 


aAUp  t2) 


3t 


j  /  /• 

-1-  -  /  M|h!|2  S*  V  Vj  [*!  -  bj  1„J  VJ  Qj1  a  -  T  bj)g  .  0 

1  0 

(3.5.1 -17a) 


3A(t1,  t2) 


3-r, 


=  f  JO){|h2!2  a*  Q2:  V2  -  b2  1M]  V*  Q-1  a  -  T  b^  =  0 

J0  “  “  ^ 

(3.5.1 -17b) 


where  for  notational  simplicity,  we  have  suppressed  the  frequency  dependency. 
The  two  likelihood  equations  can  be  further  simplified  using  the  following 
relations : 


Si/(S2  +  N)2 
1  +”61  S1/(S2  +  N) 


(3.5.1-17c) 


Sj/(S2  +  N) 

al  =  1  +  Gj.  SX/(S2  +  N) 


(3.5.1 -17d) 


tr 


-*1  <=2  tr«v2  vl)*  (V2  V*l> 


(3.5.1-17e) 
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and 


Gx  x  M 
xlx2 


X^(oj)  X^(oj) 

- r - 


Y x (tu)  Y£(<i)) 

- T - 


are  the  cross  power  spectra,  respectively,  for  X.  and 

Substituting  Equations  (3.5.1-17g  and  f)  in  Equation  (3.5.1 -17c )  yields 
the  optimum  two-target,  two-sensor  spectral  shaping  filters: 


1  + 


s1/(s2+n)2 


S2/N  |V*  V-, 
T  ~ry  N  ~2 


Sx/N 


(3.5.1-18c) 


and  similarly. 


s2/(s:+n) 

/  S1/N 

1  +  \2  '  I  +T'S"/N 


(3.5.1 -18d ) 


Note  that  the  filters  |h1 1  and  j h2 1  are  a  function  of  the  steering  vectors, 
and  that  the  target  of  interest  is  treated  as  part  of  the  signal  process  and 
remaining  targets  are  treated  as  part  of  the  noise  process. 


A  block  diagram  of  this  dual  channel  processor  is  shown  in  Figure  3-3. 
It  is  seen  that  the  dual  processor  is  tightly  coupled.  In  Chapters  4  and  5  we 
discuss  a  number  of  suboptimum  procedures  for  simplifying  the  complex  struc¬ 
ture  of  this  processor  in  light  of  improved  multitarget  time  delay  resolution. 
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Figure  3-3.  Optimum  Two-Sensor,  Two-Target  Time  Delay  Processor 
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Using  Equations  (3.4-1)  and  (3.4-5),  the  CRLB  of  the  time  delay  esti¬ 
mates  for  the  two-target,  two-sensor  case  can  be  found  as  follows: 


VAR(t • )  >  - l-T-  y-  ;  i  =  1,  2 

1  “  (i  -  mJ2)  °ii 


where  is  the  coefficient  of  mutual  dependence  given  by 


M 


12 


12  fj  1  i1/2 

(J11  J22 ' 


and  the  quantities  J. .  are  defined  by  (Equation  (3.4-5)) 

*  J 


B  r 

Vi 

k=l 


tr  - 


3ti  3tj  I 


From  the  relations 


Rk  =  Skl  Pkl  +  Sk2  Pk2  +  Nk  1 


(3.5.1- 19a ) 


(3.5.1-19b) 


(3.5.1-19C) 


(3.5. 1 -19d ) 


M 


and 


3R'1  -  ?-.i/9Pki  ~  9Gki  \  --1 

3T7-=  lhkil  * 


i  =  1,  2 


(3.5 .l-19e) 


one  obtains  after  some  straightforward  but  tedious  algebraic  manipulations: 
(See  Appendix  B) 
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where 


Gj_  =  (1  +  S2/N) 


S2/N  \ 

1  +  2  S2/N/  1-1  -2 1 


(3.5.1-22a) 


G2  *  (1  +  S^N) 


(3.5. 1 -22b ) 


and  the  quantities  Y^2  and  Y2  are  9lven  by  (Appendix  B): 


1  /  VN  \  2  Y'  A  cosn  iuA,„ 

Y1  ‘  1  '  \TT^h)  "  12 


(3.5 . l-23a) 


y12  =  cos  ul12  -  2  <  [I-  +  Gf^7(^’  +  -Nrni“+T-S27N) 


s1  s2/n‘ 


sin  ujA 


12 


(3.5.1 -23b ) 


and 


2  3 


/  S, /N  \c  ^ 

Y2  "  1  '  \  1  +  S,/N  /  £  Bn  cos 

A  n=0 


ojA 


12 


(3.5.1 -23c ) 


where  A^2  =  -  t2  is  the  time  delay  separation  between  targets  1  and  2. 

An’  Bn  are  de^inec*  ^  the  following: 


Sl/N 

1  +  s2/N 

T  +  Gj  S1/(S2  +  N) 

1  +  2  S2/N 

(3.5.1 -24a) 
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,  Sl/N  S2/N 

Xl  =  "4  1  +  G1  S1/(S2  +  N)  1  +  2  S2/N 


a2  =  1  -  Aq 


(3.5. 1 -24b ) 


(3.5.1-24c) 


A3  =  "Ai 


B0  =  4 


S2/N  ]  1  +  Sj/f 

1  +  G2  S2/(S1  +  N)  1  +  2 


(3.5. 1 -24d ) 


(3.5.1 -24e ) 


,  Sl/N 

B1  =  "4  1  +  G2  S2/(S1  +  N) 


1  +  2  Sj/N 


(3.5.1-24f ) 


b2  =  i  -  bq 


(3.5. 1 -24g ) 


B3  =  -B1  * 


(3.5. 1 -24h ) 


Note  that  if  both  targets  have  identical  power  spectra  (S^  =  S2)  and  identical 
time  delays  (t^  =  t2)  then  it  can  be  shown  that 


1  +  2  (S2/N) 

(1  +  s2/n)2 

1  +  2  Sl/N 
(1  +  Sx/N)2 


(3.5.1 -25a) 


(3.5.1-25b) 


(3.5.1 -25c ) 
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Therefore,  using  Equations  (3.5.1-25a-c)  in  (3.5.1 -20a-c ) ,  we  obtain  J,,  = 

2  2  ^  ^ 

^22  =  ^12’  Hence»  Mi2  =  ^12^11^22^  =  **  Thus>  the  time  delay  variance  is 

infinite  indicating  the  inherent  inappropriateness  of  a  two-target  formulation 

to  a  one-target  problem. 

Another  observation  is  that  letting  the  interference  power  be  zero 
(i.e.,  S£  =  0),  Equation  (3.5.1-20a)  reduces  to  Equation  (3.5.1-16a),  the  time 
delay  CRLB  for  the  single  target  case.  Note  that  because  of  the  following 
inequality: 


O  (S,/N)2 

^  T  +~ 2  Sj/N  ^  *  (3.5.1-26) 

we  conclude  that  the  time  delay  CRLB  of  the  two-target  case  (Equation 
(3.5.1-21a))  is  always  larger  than  the  time  delay  CRLB  of  a  single  target  case 
(Equation  (3.5.1-16)).  For  convenience,  let  S  =  be  the  target  power 
spectrum  and  I  =  S2  be  the  interference  power  spectrum,  then  one  can  define  a 
degradation  ratio  of  a  two-target  CRLB  over  a  single  target  CRLB  as 


(3.5.1-27) 
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Note  that  R^  _>  1  and  that  R^  is  a  function  of  SNR,  interference-to -noise  ratio 
(INR),  and  relative  time  delay  separation. 

For  the  purpose  of  illustrating  the  two-sensor,  two-target  CRLB  perform¬ 
ance  behavior,  we  assume  a  one-sided  power  spectrum  of  100  Hz  bandwidth 
centered  at  50  Hz.  Furthermore,  we  assume  that  signal,  interference,  and 
noise  processes  have  identical  bandwidths.  Figures  3-4  and  3-5  present  the 
multitarget  shaping  function  y^*  and  y^,  respectively,  as  a  function  of  fre¬ 
quency  for  a  number  of  interference-to-signal  ratios  (ISRs).  Two  effects  can 
be  seen:  (1)  the  function  y^(uj)  reduces  the  low  frequency  band  contribution 
to  the  CRLB  while  the  function  y^uj)  reduces  the  high  frequency  band  contri¬ 
bution  to  the  coefficient  of  mutual  dependence;  (2)  as  ISR  decreases,  both 
Yj(oj)  and  y^co)  approach  unity.  Figure  3-6  presents  the  two-sensor,  two- 
target  degradation  ratio  {Equation  (3.5.1-27))  as  a  function  of  time  delay 
separation  for  a  number  of  ISR  values.  Note  that  identical  signal,  inter¬ 
ference,  and  noise  power  spectra  of  100  Hz  bandwidth  are  used.  As  shown  in 
Figure  3-6,  the  degradation  ratio  is  oscillatory  as  a  function  of  time  delay 
separation,  which  is  plotted  in  terms  of  the  noise  inverse  bandwidth.  The 
oscillatory  behavior  decreases  as  separation  increases.  Two  major  peaks  are 
observed.  The  first  occurs  when  time  delay  separation  approaches  zero  and  the 
other  occurs  at  4  times  the  inverse  bandwidth.  Note  that  the  first  peak  goes 
to  infinity  at  zero  while  the  second  decreases  as  ISR  decreases.  Thus, 
targets  with  identical  spectrum  and  spatial  location  yield  the  worst  estimate 
since  they  are  not  separable  in  frequency  and/or  space.  In  order  to  see  the 
effect  of  different  signal  and  interference  spectra  on  the  resulting  esti¬ 
mate,  we  fix  the  SIR  but  vary  the  interference  spectrum.  We  choose  a  noise 
bandwidth  of  5000  Hz  (B  =  5000).  The  results  are  presented  in  Figures  3-7 
and  3-8.  Two  observations  can  be  made:  (1)  degradation  ratio  is  no  longer 
infinite  at  zero  time  delay  separation  with  non-identical  signal  and  inter¬ 
ference  spectra,  and  (2)  comparison  between  Figures  3-6  and  3-7  shows  that  the 
degradation  ratio  as  a  function  of  time  delay  separation  is  bandwidth 
independent. 


3. 5. 1.3  Case  3:  One  Target  and  M  Sensors  (J  =  1,  M  =  M).  For  M  sensors 
and  with  Qk  =  I,  the  parameter  set  consists  of  M-l  independent  time  delays. 
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Figure  3-5.  the  Function  y12(o>)  at  Different  Interference-to-Signal 
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For  convenience  we  select  the  intersensor  time  delay  as  the  minimum  time  delay 
set.  If  we  denote  this  parameter  set  by  e  =  [t^ ,  T2,  ....  then  the  opti¬ 

mum  estimate  of  this  time  delay  set  is  given  by  the  M-l  likelihood  equation: 


00 

3A  _  f  • 


|hi I2  a*  QT1  V.(*.  -  b.  1M)  V*  QT1  a  -  Tbi 


doj 

2tT 


=  0 


(3.5.1-28) 


for  i  =1,  2,  . . . ,  M-l. 


It  was  derived  in  Appendix  C  that  the  CRLB  of  time  delay  estimate  for  every 
inter-sensor  time  delay  is  identical  and  is  given  by 


VAR(Ti)  > 


B 

r 

k-l 


i'i  1  _1 


M  2  “k  mX7S 


k'  k 


i  =  1,  2,  ....  M-l 

(3.5.1 -29a) 


or  in  continuous  form 


[CO 

f  S2(o>)/N2(u>)  .2 

MT  /  r+  M  5(c3y/N(u>)  w 

0 


dui 


-1 


(3. 5.1-29b) 


Thus,  the  CRLB  of  any  two  adjacent  sensors  improves  with  increased 
M,  the  total  number  of  sensors.  Furthermore,  Appendix  C  shows  that  the 
time  delay  estimate  between  any  two  sensors  also  has  the  same  CRLB.  In 
addition,  the  covariance  of  any  two  time  delay  estimates  is 


COV(t1,tj.) 


if  | i-j |  =  1 
if  | i-j  |  >  1 


(3.5.1-30a) 


In  other  words,  the  correlation  coefficient  is 
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* 


* 

■ '  . 

:  . 

1 


1 


4 


ij 


cov(y  t.) 


VAR(  t.  )  VAR(x,f2 

1  J 


-  \  ;  if  I i  -  J  I  =  1 
0  ;  if  j i  -  j |  >  1 


(3-5.1 -30b) 


Thus  far,  we  have  studied  only  the  optimal  time  delay  processing  of  a 
multisensor  array  system.  The  objective  of  a  passive  sonar  system,  however, 
is  to  localize  an  acoustic  source  of  interest  (i.e.,  obtaining  its  range  and 
bearing).  For  passive  localization,  a  minimum  of  three  sensors  (assuming 
omni-directional  response)  is  needed.  However,  for  a  far  field  assumption, 
two  sensors  can  yield  good  bearing  information.  In  the  following  section,  we 
study  in  some  detail  the  optimum  processor  structure  for  range  and  bearing 
estimations  from  a  three-sensor  array. 


3.5.2  Localization  Parameter  Estimation 


In  this  section  we  examine  in  detail  the  optimum  processor  structure  for 
estimating  localization  parameters  (i.e.,  range  and  bearing).  We  compare  the 
performance  of  the  optimum  processor  to  the  conventional  approach,  where 
range  and  bearing  are  obtained  from  the  measurement  of  two  time  delay  pairs. 

We  discuss  briefly  the  optimality  of  the  direct  range  and  bearing  estimation 
approach  and  of  the  indirect  time  delay  estimation  approach. 

Consider  a  general  three-sensor  array  geometry  shown  in  Figure  3-9. 

We  are  interested  in  obtaining  an  optimum  estimation  of  range  and  bearing 
using  outputs  from  a  three-sensor  array.  From  Appendix  D,  we  note  that 
the  range  and  bearing  information  is  contained  in  the  two  incremental 
time  delays.  Let  denote  the  time  delay  between  sensor  1  and  2,  and 
let  X2  denote  the  time  delay  between  sensors  2  and  3.  Then  one  can  write 

r  -  %  r^  +  L?  -  2rL,  cos(tt  -  e) 

T1  =  °2  -  D1 - c -  (3.5.2-1) 


fl 
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(3. 5. 2-2) 


Jr^  +  l~  -  2rL9  cos(e  -  <J>)  -  r 

t2  =  Ds  -  °2  =  1 ^ - 

where  Di ;  i  =  1,  2,  and  3  are  the  propagation  delays.  The  steering  matrix 
is 

t  -JXti  J<V2  ) 

Vk  =  diag  j  e  *  1  1  e  *  L  j  . 

Therefore,  Equations  (3.5.1-10a)  and  (3.5.1-10b)  become 

Pk  "  vk  *M  \ 

3Pk  .  w  3$_  v* 

3ei  J“k  Vk  3ei  Vi< 

where  el-  is  any  time-delay  related  parameter  of  interest,  and  the  matrix 
$  is  given  by  (see  Equation  3.5.1-8d) 


(3. 5. 2-3) 

(3.5.2-4a) 

(3.5.2-4b) 


y 


Figure  3-9.  A  Three-Sensor  Array  System 
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$  = 


-Tl  -(T!  +  T?) 

r2 


rl  "V  ll  T  l2; 

0  -Tr 


(Ti  +  T?) 


1  T  2'  2 


0 


(3.5.2-4C) 


The  likelihood  equations  for  range  and  bearing  can  be  obtained  from  Equation 
(3.5.1-lla)  as 


B 


k=l 


l\i2  <S  Q4  \  (f  -  iM)  Qk4  a,  -  Sr 


(3.5.2-5a) 


and 


3A(r,e)  .. 


3e 


■  E  n  'ski2  <*  54  vk  (||  -  b'9 1„)  vk.  q-1  s, .  b'( 

k=l 


=  0 


(3.5.2-5b) 


where  from  Equation  (3.5.1 -lib)  we  have 


b'r  ■  >k  \  W  Vf) 


(3.5.2.5c) 


b0  =  ak 


«(<>*  vk  S  ve)- 


(3.5.2 -5d) 


For  a  single  target  in  uncorrelated  noise,  it  is  easily  verified  that  the 
biases  are  zero,  i.e.,  and  b  =  b  =  0.  Furthermore,  for  convenience  we  write 

-  ~  i  r  9 

lhkl  =  l\l  and  ^k  =  1,e*»  assuming  spectrally  identical  but  spatially 

uncorrelated  noise  power  spectrum  for  each  sensor.  Therefore,  the  range  and 
bearing  likelihood  equations  reduce  to 

■  E  J"k  |hk|2  Vk  If  Vf  2k  -  0  (3.5.2.6a) 

k=l 
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(3.5. 2 -6b ) 


where  from  Equation  (3.5.1-4b)  one  obtains 
?  Sk/Nk 

lhkl  =  1  +  3  $k/Nk  *  (3.5.2-6c 

Using  the  definitions  of  Vk,  $  and  the  resulting  derivatives  9<J>/9r  and  9$/9e, 
the  likelihood  equations  can  be  manipulated  to  yield  (Appendix  F) 


00 

=Jf| r  f  |hU)|2  G(0);  t2)  do,  =  0  (3.5.2-7a) 

-00 


*  hh  f  ih(“>i2  G(“:  Ti,  t2>  d“  * 0 


where 


(3.5. 2-7b) 


G  (o>; 


Tl. 


t2) 


JUJTi  JCOTo 

=  G2^(o))e  +  G7?(oj)e 


32' 


+  G31(ai)e 


jw(t 


1+t2) 


(3.5.2-7c) 


and  G^u)  =  Ta^a*  is  the  estimated  cross-power  spectrum  between  sensor  i 
and  j.  Now  we  define  the  joint  parameter  ambiguity  function  by 


R  (t^  t2)  = 


/  |h(u,)|2  G((d;  t1s  t2)  % 


=  J  I h(co)  | 2  G21(uj)e  1  +  G32(io)e  2  +  G31(oo)e 


jd)(T1+T2) 


du) 

7T 
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=  ^21 ^  T1 ^  +  ^32  ^  t2  ^  +  ^31 ( Tl  +  t2^  ’ 


(3. 5. 2-8) 


Then  the  jo'nt  estimation  of  range  and  bearing  can  be  realized  by  an  open  loop 
processor  described  in  Figure  3-10.  In  the  literature,  this  is  known  as  a 
focused  beamformer  because  the  optimum  estimate  of  range  and  bearing  is 
obtained  when  the  beam  is  focused  on  the  target. 

Note  that  the  optimum  estimate  of  range  and  bearing  (r*,  e*)  must  cor¬ 
respond  to  the  time  delay  pair  (t£,  t£)  which  defines  the  peak  of  the  joint 
time  delay  ambiguity  function  R(x^  Xg).  Therefore,  an  equivalent  theoretical 
approach  is  to  first  seek  the  optimum  time  delay  pair  (if,  x£)  and  then  trans¬ 
form  it  to  the  corresponding  range  and  bearing.  For  practical  applications, 
however,  it  is  sometimes  more  convenient  and  simpler  to  search  and  track  in 
the  time  delay  parameter  than  the  range  and  bearing  parameters  since  the  ambi¬ 
guity  function  is  symmetrical  and  uni-modal  (for  high  SNR)  in  the  time  delay 
variables.  More  importantly,  for  a  low  SNR  environment  where  a  T-second 
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Figure  3-10.  Optimum  Range  and  Bearing  Estimation  from  Three-Sensor  Arrays  (Focused  Beamformer) 
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observation  MLE  processor  fails  to  be  efficient,  the  symmetric  property  of  the 
time  delay  ambiguity  function  allows  a  simple  tracking  filter  design  which 
could  provide  optimum  estimates  by  increasing  the  effective  coherent  integra¬ 
tion  time.  Unfortunately,  the  resulting  non-linear  transformation  from  time 
delays  to  range  and  bearing  for  the  practical  approach  renders  the  otherwise 
Gaussian  noise  process  to  be  non-Gaussi an. 

Note  the  focused  beamformer  implementation  shown  in  Figure  3-10  requires 
a  two-dimensional  (2-D)  peak  detector.  However,  a  time  delay  approach  using 
only  two  1-D  null  detectors  can  be  realized  without  any  loss  of  performance. 
This  can  be  seen  as  follows.  Rewriting  Equations  (3.5.2-9a  and  b)  in  terms  of 
the  time  delay  variables  and  using  the  chain  rule  of  differentiation  yields 


3A(r,e) 

3r 


3A(r,e) 

3e 


=  T 


=  T 


'3R(t,  To  ) 

3T1  + 

3R(tu  t2) 

3t2 

3T1 

~w + 

3t2 

3r 

'SRUj  t2) 

SR^  t2) 

3t2  1 

3e 

3e 

=  0 


*  0  . 


(3.5.2-10a) 

(3.5.2-10b) 


Using  the  simple  relations 


2  2 

L,  Lf  sin  c  9 

Ti  =  -  —  cos  e - ^ -  and 

2  2 

l_2  Lo  sin  e 

t2  =  -  -  cos  e  +  —255 - 


as  shown  in  Equations  { D-2a )  and  (D-2b),  respectively  in  Appendix  D,  we 
observe  that  the  derivatives  of  and  t2  w.r.t.  range  are  zero  only  for  infi¬ 
nite  range  and  the  derivatives  w.r.t.  bearing  are  zero  only  at  the  endfire 
direction.  Therefore,  for  all  other  situations,  the  necessary  and  sufficient 
conditions  for  Equations  (3.5.2-10a)  and  (3.5.2-10b)  to  be  true  are 
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3R(t1s  t2) 


3t, 


=  0 


(3.5 .2 -10c ) 


3R(ti  t2) 


3t, 


=  0  . 


(3.5.2-10d) 


Thus,  using  Equation  (3. 5. 2-8)  in  (3.5.2-10a-d)  yields 


3R(T1,  t2}  =  T  _3_ 
3t,  ~  3t, 


R21 ( T1 )  +  R31^t1  +  t2^ 


T  a  ?  ,  f  jw(T,+T?)  I 

=  Tir  "3r""  /  MWH  ^(wje  +  G^(u))e  I  doo 

•  00 

(3.5.2-10e) 


and 


3R<T1.  T2>  T  3 

3  T0  "  3  To 


R32(t2}  +  R31(t1  +  t2} 


jlof^+Tg) 


=  7?  "3t~  J  M^l2  ®32 2  +  x  J  du» 

^  -oo 

(3.5.2-10f ) 

where  we  have  used  the  fact  that 


3R32^t2^  3R21^t1^ 

3r^i  3tI 


=  0  . 


(3.5.2-10g) 


An  optimum  three-sensor  array  processor  using  two  correlators  is  shown 
in  Figure  3-11. 

We  next  calculate  the  CRLB  for  the  optimum  three-sensor  array  processor. 
First  we  calculate  the  CRLB  for  the  optimum  time  delay  estimate  and  compare  it 
to  the  conventional  approach.  Then  the  CRLB  for  the  localization  parameters 
is  obtained  by  relating  it  to  the  time  delay  estimate. 
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Figure  3-11.  Optimum  Range  and  Bearing  Estimation  from  Three-Sensor  Arrays  Using  Two  Correlators 
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From  Appendix  C,  we  obtain  the  elements  of  the  Fisher  Information  matrix 
for  optimum  one-target  three-sensor  time  delay  estimation  as: 


Jij  * 


B  / 
/  3$  3$  \  V'  2  I 


S^/N^ 

VNk 


3  Sk/Nk 


(3.5.2-11) 


From  Equation  (3. 5.2 -4c )  we  find 


0  -1  -1 

M_  =  i  o  o  •  — 

3ti  U  *  3t5 

1  1  0  0  c 


0  0-1 

0  0-1 

1  1  0 


(3.5.2 -12a ) 


and  therefore 


/  3$  3$_\ 

l  3xl  3t1  / 


(3. 5.2 -12b ) 


/  3<I>  3$_  \  _  , 


(3.5.2-12C) 


-  (%%)•■*• 

From  Appendix  B,  the  CRIB  for  a  two-parameter  joint  estimation  is: 


(3.5.2 -12d ) 


VAR(^)  >  - X-  X 

1  (i  -  m;2)  Jn 


(3.5.2-13a) 


VAR(t-)  >  - p — 


(1  -  MJ2)  °22 


(3. 5.2 - 13b ) 
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where  M^»  t,ie  coefficient  of  mutual  dependence,  is  given  by 


n 


J12 


12  "  (Jn  'J22)Js  ’ 


(3-5.2 -13c ) 


Now  using  Equations  (3.5.2-11)  and  (3.5.2-12a-d)  or  Equation  (3.5.1-29a),  we 
obtai n 


VAR(?t)  =  VAR(t2)  > 


B 


3S“k 

k=l 


2/  ^k 


1  +  3  Sk  Nk> 


-1 


(3.5.2-14) 


Note  that  for  the  general  case  of  M  sensors  it  was  shown  in  Appendix  C  that 
the  CRLB  of  the  incremental  time  delay  estimates  is  given  by 


VAR(f.)  > 


B 


2 1  Sk/Nk 


L  k 


L/^kyi  +  m  sk/Nk/J 


-l 


i  =1,  2,  ...,  M-l  . 


(3.5.2-15) 


Thus,  the  time  delay  estimate  improves  with  the  increased  number  of  spatial 
sensor  observations. 


N 


F« 


We  now  discuss  the  conventional  approach  to  estimating  range  and  bear¬ 
ing.  A  conventional  three-sensor  array  signal  processor  is  shown  in  Fig¬ 
ure  3-12.  In  the  conventional  approach,  only  two  cross-power  spectra  are 
processed.  For  this  case,  the  CRLB  of  the  time  delay  estimate  is  given  by 
(see  Equation  ( 3. 5 .1 -16a) ) 

B 


VARC(?1)  =  VARc(t2)  > 


to,. 


sk'Nk 


L  k 


Tl  *  U  +  *  Sk/Nk 


-1 


(3.5.2 -16a ) 
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Therefore,  the  time  delay  variance  ratio  of  the  optimum  to  the  conventional  is 


(3.5.2-165) 


and  for  a  flat  signal  and  noise  power  spectra.  Equation  (3.5. 2 -16b )  is  simpli¬ 
fied  to 


VAR(^) 

VAR-tfJ  = 


1  +  3  S/N  \ 
ITT  S/N  /  ‘ 


(3.5.2-16c) 


Thus,  at  low  SNR,  the  CRLB  of  the  optimum  processor  is  improved  by  a  factor  of 
2/3  (1.8  dB)  with  respect  to  the  conventional  approach. 


Because  the  two  time  delays  are  correlated  due  to  common  noise  channels,  a 
more  meaningful  approach  is  to  compare  the  one  sigma  error  ellipse  area 
between  the  optimum  and  the  conventional.  From  Appendices  C  and  E,  we  obtain 
the  two  time  delay  covariance  matrices  as 


P 


o 


1 


-1 

~7 


-1 

T 


VAR(t1) 


(3.5.2-16d) 


P 

c 


1 

-S/N 

1  +  2  S/N 


-S/N 

1  +  2  S/N 
1 


VARC(^) 


(3.5.2-16e) 
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where  PQ  and  Pc  denote,  respectively,  the  optimum  and  conventional  covariance 
matrices.  Therefore,  the  ratio  of  their  area  is 


-|P0I  2 


C  IT  IP, 


(1+2  S/NT 
+  3  S/N)(l  +  s/n; 


VAR(^) 

varc(t1; 


(3. 5.2 -16f ) 


Using  Equations  (3.5.2-16d-e)  in  Equation  (3.5.2 -16f )  yields  the  desired  result 


1 1  1  +  3  S/N 
I  "1  +  S/N 


(3. 5 . 2 -16g ) 


Thus  Equation  (3.5.2 -16g )  states  that  in  a  low  SNR  environment  (S/N  «  1), 
the  error  ellipse  of  the  optimum  approach  has  an  area  equal  approximately 
to  one  half  the  conventional  approach.  This  represents  a  substantial  loss 
in  performance  by  the  conventional  approach.  On  the  othe^  hand,  in  a  high 
SNR  environment  (S/N  »  1)  the  conventional  processor  approaches  the  optimum. 

The  CRLB  on  the  localization  parameters  for  the  focused  beamformer  can  be 
obtained  as  follows.  Defining  the  Fisher  Information  matrix  by 


F  F 
rr  re 


Fa„  Fee 

er 


then  the  CRLB  for  the  range  and  bearing  estimates  is: 


VAR(r) 


(1  -  Mj2)  rr 


(3.5.2-17a) 
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4 


VAR(e)  = 


1 


(1  - 


M12> 


1 

Fee 


where 


re 


12 


<Frr  Fee)' 


Again  from  Appendix  C,  we  find: 

B 

P  f  /  c$  30  \  V'  2 

Frr  "tr  (  3r  dr)  2-j  u k 

k=l 


/  a« 

a$  \ 

B 

2  / 

4»l  \ 

\  ar 

3e  / 

k=l 

“k  \1 

1  +  3  VNk  / 

Sk/Nk 

1  +  3 . S.7N7 


and 


VNk 


1  +  3  Sk/N 
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(3.5.2 -17b ) 


(3.5.2 -17c ) 


(3.5. 2 -18a ) 


(3.5. 2 -18b ) 


(3.5.2-18c) 


(3.5.2-19a) 
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and  using  Equations  (3.5.2-12a-d),  we  obtain 


/3£3£\_  4/!ll  +  !!l  ^2  ^2  ^2  ^2  )  (35219b) 

Ur  3 e/  4  \  3r  3e  3r  3e  3r  3e  3r  3©  /  (3. 5. 2-19b ) 


and  similarly 


/3$3$\_  .  /^l\  2  ^1  +  (  8t2  \  2 

\  3r  3r/  “  \  3r  /  3r  3r  \  3r  / 


(3.5.2-19c) 


tr/3ii£\  ,  .4[(!li)  2  +  !li!l2  +(!!i)2|  (3.5.2-19d) 

\3e  3e/  4  [  \  3e  /  3e  3e  \  3©  /  * 


Using  Equations  (3. 5. 2-1)  and  (3. 5.2-2)  and  evaluating  the  derivatives 
at  the  true  range  and  bearing  (R,  B),  we  obtain  the  following: 


l  fl2  si n(B  -  <J>) 
7  \ - R 


( 3  -  5 . 2 -20a ) 


(3.5.2 -20b ) 


(3.5.2 -20c ) 


^®|r,b  c 


(3.5.2 -20d ) 
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For  simplicity  we  assume  a  special  case  where  the  sensor  arrays  are  co-linear 
with  equal  separation  (L  =  =  l_2)  and  the  target  is  at  broadside.  Thus,  sub 
stituting  Equations  (3.5.2-20a-d)  in  Equations  (3.5.2-9b-d)  yields: 


(3.5.2-21a) 

(3-5.2 -21b) 


and 


(3.5. 2 -21c ) 


Therefore,  the  CRLB  of  range  and  bearing  evaluated  at  the  true  range  and  bear 
ing  is  given  by: 


VAR(r)  >  - - — 5 — 

~  (i  -  mJ2) 


-tr 


w-1 


(  3r)  |  S  “k  (l  +  3  Sk/Nk  /j 


(3.5.2-22a) 


VAR(e)  > 


(1  -  M12)  I 


1  .  \(*)*}  V  z(  <><  VI 
|  ”tr  [  (  3©  /  J  k  v  1  +  3  Sk/Nk/| 


-1 


(3.5.2 -22b) 


And  because  of  Equation  (3.5. 2-21  a) ,  range  and  bearing  estimates  are 
uncorrelated.  Using  Equation  (3.5.2-14)  these  can  be  expressed  in  terms  of 
the  variance  of  the  time  delay  estimate  as  follows: 
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VAR(r)  > 


3  VAR(t1) 


(3.5.2-23a) 


VAR(e)  > 


3  VAR  (Tl) 


(3.5. 2 -23b ) 


For  the  special  case  of  a  co-linear  array  and  a  broadside  target,  we  find 


VAR( r )  >  3c2  (£)4  VAR(tj) 


(3.5.2-24a) 


VAR(e)  >\(if  VAR(t1)  . 


(3.5. 2 -24b ) 


Taking  the  ratio  of  Equation  (3.5.2-24a)  to  Equation  (3.5.2-24b)  yields 


VAR(r)  >  12  K  VAR(e)  , 

l/ 


(3.5. 2 -24c ) 


which  relates  the  range  variance  to  the  bearing  variance.  Thus,  we  obtain  the 
well  known  results  that  the  variance  of  the  range  estimate  is  proportional  to 
the  fourth  power  ratio  of  the  true  range  to  base  line  length,  and  the  variance 
of  the  bearing  estimate  is  inversely  proportional  to  the  square  of  the  base 
length. 

Note  that  the  range  and  bearing  variances  of  the  focused  beamformer 
approach  (Equations  (3.5.2-24a-b))  agree  with  Equations  (G-15)  and  (G-16)  of 
Appendix  G,  which  were  obtained  via  a  geometric  mapping  from  time  delay 
measurements.  Therefore,  the  one-step  focused  beamformer  approach  and  the 
two-step  time  delay  approach  yield  identical  statistical  performance. 


Let  SQ  and  SQ  denote  the  area  of  the  range/bearing  one-sigma  localization 
error  ellipse  for  the  optimum  and  the  conventional  approach,  respectively. 

Then  from  Equations  (3.5.2-24a)  and  (3.5.2-24b)  one  obtains 
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(3.5. 2-25a) 


SQ  =  TrVvAR(r)  VAR(e) 

On  the  other  hand,  for  the  conventional  approach,  it  can  be  shown  (Appendix  G) 
that 


■  W 


2c2\Vu~ 


+  3  S/N)(l 


+  S/N) 

VW) 


VARc(t1) 


(3.5. 2-25b) 


Therefore,  using  the  relation  given  in  Equation  (3.2.5-l6d),  the  ratio 
of  the  error  ellipse  area  is 


(3.5.2-25 c) 


which  has  the  same  ratio  as  in  time  delay  estimation  (see  Equation  3.5.2-16g). 
Thus  Equation  (3.5.2-25c)  implies  that  the  optimum  processor  yields  a  one- 
sigma  localization  error  ellipse  which  is  approximately  one  half  (1//T)  that 
of  the  conventional  approach  in  a  low  SNR  environment.  This  improvement  comes 
directly  from  a  better  bearing  estimation  using  the  optimum  approach.*  Note, 
it  can  be  shown  that  the  ranging  performance  is  identical  between  the  conven¬ 
tional  and  the  optimum  approach.  (See  Appendix  G.) 


*It  was  pointed  out  by  Dr.  0.  Iannielo  of  the  Naval  Underwater  Systems  that 
optimum  range  and  bearing  estimation  can  also  be  achieved  using  the  conven¬ 
tional  approach.  However,  one  must  provide  separate  and  different  spectral 
shaping  filters  for  range  and  bearing  estimation. 
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3.5.3  Power  Spectral  Estimation 


Our  discussion  thus  far  assumes  that  all  target  power  spectra  are  known. 
This  is  one  of  the  strongest  assumptions  we  have  made  in  studying  the  optimum 
signal  processor.  The  resulting  processors  contain  spectral  shaping  filters 
which  are  functions  of  the  known  target  power  spectra.  In  an  actual  implemen¬ 
tation,  the  power  spectra  must  be  either  known  a  priori  or  estimated.  In  this 
section  we  briefly  examine  the  methodology  of  spectral  estimation  in  a  multi¬ 
sensor,  multitarget  environment  and  the  relations  between  power  spectral 
estimation  and  time  delay  estimation. 


Thus,  we  seek  the  estimate  S^.,  the  signal  power  spectral  level  of 
target  j  at  frequency  k.  Recall  from  Equation  (3.2-6)  that  the  spectral 
levels  are  contained  in  the  observation  covariance  matrix  R^.  Let  S  be  a 
column  vector  of  all  the  unknown  spectral  levels.  Thus, 

S  =  (S^  S21  .  .  .  SB1;  S12  S22  .  .  .  Sg2;  .  .  .  ;  S2J  .  .  .  SBJ)  is  a 

JB  dimension  vector.  Therefore,  from  Equation  (3.2-9)  the  likelihood  equa¬ 


tion  for  Skj. 


is 


3A(S) 


3Ak(S) 


3S 


kj 


=  0; 


j  =  1,  2, 
k  =  1,  2, 


•  ,  J 

•  ,  8 


(3. 5. 3-1) 


Note  that  the  likelihood  equations  are  decoupled  in  frequency.  This  implies 
that  each  equation  can  be  solved  separately.  Now  from  Equation  (3.2-10) 
we  obtain 

S)  /  i  3Rk  \  *  SR:1 

—  =tr(Rk  9S^j+ -^k  3S^7  (3. 5. 3-2) 

where  Rk  and  Rk*  are  given  by  Equations  (3.5.1-2a)  and  (3. 5. 1-3).  Hence, 
the  derivative  of  Rk  and  Rj^  w.r.t.  Skj-  is  "iven  by 


3Ak( 

~^ST 
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Substituting  Equations  (3.5. 3-3b)  and  (3. 5. 3-4)  into  (3. 5. 3-2),  one  obtains 


3A(  S) 


«kUi2 


1  *  Gkj  Skj/Nk3  (1  +  Gko  W  Nkj 


(3. 5. 3-5) 


Solving  Equation  (3. 5. 3-5)  yields 


;  l\i  GkU'2  \±  . 


k  =  1,  2,  . 

j  -  i,  2,  . 


.  ,  B 
•  ,  J 


(3. 5. 3-6) 


Thus,  in  order  to  estimate  the  J  target  power  spectra,  a  total  of  JB  equations 

★ 

must  be  solved  if  the  optimum  steering  vector  is  k  nown.  If  it  is  not 

known,  then  the  J(M  -  1)  time  delay  equations  must  be  solved  simultaneously 
with  the  JB  spectral  equations.  However,  in  practice  the  target  spectrum  can 
be  modeled  with  a  considerably  smaller  set  of  unknown  parameters.  Therefore, 
sampling  the  spectrum  at  appropriate  frequencies  should  provide  sufficient 
information  to  estimate  the  unknown  spectral  parameters. 

Note  that  Equation  (3. 5. 3-6)  is  an  unbiased  spectral  estimator.  This  can 
be  shown  easily  as  follows. 

Taking  the  expected  value  of  Equation  (3. 5. 3-6)  yields 


^  -  £ 
xj  KJ 

Gk'j  E<^2k*>  Gk1 


(3. 5. 3-7) 


TR  67  7 


But  from  Equations  (2.3-3a)  and  ( 3 . 2 -7a) : 

E<2k  2k*)  -  SkJ  Vy  4*3  *  SkJ  Qkj  • 

Substituting  Equation  (3. 5.3-8)  into  (3. 5.3-7)  yields 


E‘skj) 


s  .  Nk.i  ^  ii 

kJ  ^ 


(3. 5. 3-8) 


(3. 5. 3-9) 


Next  we  briefly  examine  the  optimum  spectral  estimator  for  the  one-target 
and  two -target  cases. 


3. 5. 3.1  Case  1:  One  Target,  M  Sensors  (J=l,  M=M) .  For  a  single  target, 
we  have  Nkj.  =  Nk,  Qkj  =  Qk,  and  Gkj.  =  Gk  =  Qk*  Equation  (3. 5. 3-6)  becomes 


A 


(3.5.3-10) 


Furthermore,  by  letting  Qk  =  I  (i.e.,  noises  are  uncorrelated  from  sensor 
to  sensor),  then  Equation  (3.5.3-10)  can  be  simplified  to 


a  i£ski2  Nk 

k  m2  -  M 


(3.5.3-11) 


Figure  3-13  shows  the  optimum  one-target  multisensor  spectral  estimator. 


3. 5. 3. 2  Case  2:  Two  Targets,  M  Sensors  (J=2,  M=M) .  For  the  two-target 
case,  Equation  (3. 5.3-6)  becomes 
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kl 


j4l  ^1%!'  Sk2  *  \ 

r2 
Gkl 


.2  P 


3kl 


(3.5.2-12a) 


k2 


14*2  Qkla/  Skl  +  Nk 


k2 


k2 


(3.5. 2-12b) 


Recall  that 


P-1 


Qkl  *  <Sk2  *  Nk>  <Sk2  42  42  +  \  V'1 


»  +  Sk2/Nk) 


■1 


Sk2/Nk 


Qk1  Ir,  h,  Q,'1 


1  + "  v\ Vk  ^2  ^2  "k 


Simil  arly, 


(3.5.3-13a) 


qki  -  a  *  skl/Nkj 


hi** 


<k  1  *  M  Skl/Nk  5k  ^kl  4l  Qk 


(3. 5. 3-13b) 


Also 


4i  *  4*i  4i 


4*i  ^k1  4i  -  — <2'"k  '»*  "-1”  -2 


^k2^*k 

: — ^ —  (Vfci  q’1  v,,)1 

1  +  M  Sk2/Nk  k 


(3.5.3  -13c) 
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d  +  Ski/Nk)  4%  Q-1  42  -  $klf k  (4*2  Qk1  42)2 

1  *  N  $kl/Nk 


(3.5. 3 -13d ) 


*  -1  *  1 

For  simplicity  assume  Qk  =  I,  then  hi  V  4l  '  Vk2  V  h2  =  "■  Now 
the  following  relations  can  be  obtained: 


"  n  *  sk2/Nk) 

1  *  M  Sk2/Nk 

”  n  +  ~su^k> 

1  *  N  Skl/Nk 


14*1  q,1!  a,!2  >  (i  *  ?k2/Nk)2  4*1  a. 


Sk2^k 

1  +  M  ^k 


(3.5.3-14a) 


(3.5. 3-14b) 


(iiy  <&*> 


(3.5.3-146 ) 


S  /N  ; 

14*2  Qk1  Ski2  =  (1  *  4i/Nk)2  4*2  2k  -  -  kl> k-  (4*1 42>  til  Sk>  ' 

i  ♦  n  ^kl/Nk 

(3.5. 3-14d ) 


Substituting  the  above  into  Equation  ( 3.5.3-12a-b)  and  simplifying  yields 
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Skl  =  Ak2 


Ak2  l^kl  Sk 


Bk2  hz  ^kl2 


(3.5. 3-15a ) 


Sk2  =  Akl 


Akl  ^k  '  Bkl  ^kl  ^k  1 2 


(3.5. 3-15b) 


where  A^.  and  are  defined  by 


1  +  MSki/Nk 


ki 


i  -  1,  2 


(3.5. 3-15c ) 


and 


B, 


^i^k 


ki  ’  1  +  M  Ski^k 


V  V 

4l  \2' 


i  =  1,  2  . 


(3.5. 3-15d ) 


The  optimum  two-target  M  sensor  spectral  estimator  is  shown  in 
Figure  3-14. 

We  have  already  shown  that  the  optimum  spectral  estimator  is  an  unbiased 
estimator.  In  the  remainder  of  this  section  we  briefly  discuss  the  estimator 
performance  bound.  There  are  two  cases  of  particular  interest:  (1)  power 
spectral  estimation  with  known  time  delay,  and  (2)  joint  time  delay  and 
spectral  estimation. 

3. 5. 3. 3  Case  1:  Power  Spectral  Estimation  with  Known  Time  Delay. 
Consider  the  two-target  M  sensor  problem.  The  optimum  power  spectral  estima¬ 
tor  is  given  in  Figure  3-14.  The  spectral  CRLB  can  be  derived  as  follows. 
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For  each  frequency  co^,  we  have  the  relation 


VAR<Skj>  ■  'F %  ■  J  *  2 


(3.5.3-16) 


where  F  is  the  2x2  Fisher  information  matrix  whose  j£,  element  is  given 
by 


a2A(S) 
3Skj  3Sk£ 


/  _i  -1  3^k 

strVk  ^7 Rk 


*  -tr 


/  3^3^  \ 
\  3Skj  3Skl  / 


(3.5.3-17) 


where  Equations  (3. 5. 3-2)  and  (3.4-4)  have  been  used.  Using  Equations  (3.5.3-3a) 
and  (3.5.3-3b),  Equation  (3.5.3-17)  reduces  to 


F  -  tr(  Gkj  ^3  W"kj  ,J  v* 

J1  \  f1  *  Gkj  vv2 


M  \i'2/iik.i 
<!  *  Gkj  V"kj> 


(3.5.3-18a) 
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Note  for  j  =  i,,  Equation  (3.5. 3-iSa)  further  reduces  to 


F  .  VA.1 

L1  +  Gkj  skj/Nkj 


(3.5. 3 -18b ) 


Thus,  the  spectral  CRLB  is 


VAR(S.  .)  = 


(1  -  Mto)  JJ 


(1  -  ) 


(s  • 


j  -  1,  2 


(3.5. 3-19a) 


where 


11  r22'  2 


l"k*2  ft  «ll 


Jkl  k2' 


(3.5. 3-19b) 


is  the  spectral  coefficient  of  mutual  dependence.  Note  that  for  a  single  tar¬ 
get  case,  the  spectral  CRLB  is 
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VAR(Sk)  =  - 


-1 


+ 


2 


(3.5.3-20) 


where  Gk  =  V£  is  known  as  the  array  gain. 

3. 5. 3. 4  Case  2:  Joint  Time  Delay  and  Spectral  Estimation.  For  sim¬ 
plicity,  we  assume  a  single  target  and  two  sensors.  Therefore,  the  unknowns 
are  t,  the  time  delay,  and  $k,  the  spectral  level  at  frequency  oj^.  Note  that 
because  the  spectral  likelihood  equation  is  independent  for  each  fre¬ 
quency  wk,  we  only  need  to  consider  the  joint  estimate  between  t  and  Sk.  Now 
let  e1  =  t,  and  e2  =  Sk’  then  the  J°int  time  delay  spectral  CRLB  evaluated  at 
the  true  parameter  values  is 


VAR(e.)  >  (F"1)ii  ,  1  •  1,  2 


(3.5. 3-21  a) 


where  the  ij  element  of  the  Fisher  information  matrix  is  given  by 


(3.5.3-21b) 
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/  VNk  \2 
V  +  Gk  VNk  / 


(3.5. 3-21d ) 


Finally,  the  cross  term  is 
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JdZh(S,T)\ 
F12  =  "E  \  3t  ask  / 

■~(£s) 


’  J“klhkl  tr<vk  vt  vk  W 


(3.5.3-21e) 


where  relations  in  Equations  (3.4-2)  and  (3.4. -5b)  have  been  used.  But  from 
Equation  (3.5.1-13e),  we  obtain 


trl\  *1  V£  Vk  !M  Vk>  *  tr 


0  -1  [1  o 


0  e‘JU5T  1  OJ  0  e‘JWT 


0  eJWT  1  1 


'-1  -e”J(i)T 
tr  ejuT  l 


*  0  . 


(3.5. 3-21f ) 


Thus,  =  0,  indicating  that  spectral  estimates  and  time  delay  estimate  are 
uncorrelated.  Hence,  using  Equation  (3.5.3-21a)  the  CRLB  is  given  by 
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VAR(^)  -  (2^  4  |hk|2  Sk)  -1 


(3.5.3.22a) 


VAR(Sk)  = 


•R) 


(3.5.3-22b) 


which  is  identical  to  the  time  delay  estimate  with  a  known  spectrum  and  to  the 
spectral  estimate  with  a  known  time  delay,  respectively.  Thus,  we  conclude 
that  joint  time  delay  spectral  estimation  does  not  degrade  the  time  delay 
estimates  nor  the  spectral  estimates. 

It  is  interesting  and  revealing  to  show  the  explicit  dependence  of  the 
time  delay  spectral  performance  on  the  observation  time  T.  Let  S(co)  and  N(u) 
be  the  true  continuous  signal  and  noise  power  spectra.  Now  recall  from 
Equations  (2-13d)  and  (3.5.1-13d)  that  S.,  =  and 

«  n  l\  Is  l\ 

l\l  =  Sk^k^  +  2  Sk>/lV’  where  \  is  the  fr°urier  coefficient  of  the 
signal  waveform  from  T  seconds  of  observation.  Also  recall  the  relation  that 
for  sufficiently  large  T,  we  have  TSk  =  TNk  =  N(tuk),  and 

|hk|2  =  T|h(uk))2.  Consequently,  Equations  (3.5.1-22a)  and  (3.5.1 -22b)  can 
be  manipulated  as  follows  to  yield 


VAR(t)  =  I  co2  |h(wk)j2  S(u>k)  h 

T 


=  -  f 

T 


S2(cj) /N2  (cu) 

+  c  S(oj)/N(oj) 


(3.5.3-23a) 
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VAR(Sk)  -  ij 


(TNk) 

(TSk>  *  -sr 


S(CJk )  + 


(3.5. 3 -23b ) 


where  S^)  denotes  S(co)  evaluated  at  frequency  w  =  2irk/T. 

Thus,  Equations  (3.5. 3-23 a)  and  (3.5. 3-23b )  indicate  that  while  the  time 
delay  variance  is  inversely  proportional  to  T,  the  spectral  variance  is 
inversely  portional  to  the  square  of  T.  In  other  words,  the  spectral  variance 
is  more  effectively  reduced  by  increasing  the  observation  time  than  the  time 
del^  variance. 
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CHAPTER  4 

SUBOPTIMUM  REALIZATION  OF  MULTI  SENSOR,  MULTITARGET  TIME  DELAY  PROCESSOR 

4.1  INTRODUCTION 

In  Chapter  3  we  derived  the  optimum  (MLE)  multisensor,  multitarget  time 
delay  estimator.  The  result  yields  a  highly  coupled  multi-channel  processor. 
For  practical  applications,  it  is  desired  to  seek  suboptimum  realizations 
which  can  substantially  simplify  the  required  implementation.  In  this  sec¬ 
tion  we  examine  the  suboptimum  processor  based  on  a  weak  signal  in  noise 
assumption.  This  is  the  case  of  considerable  interest  since  in  passive  signal 
processing,  a  weak  signal  in  noise  represents  the  usual  environment  at  which  a 
signal  processor  must  operate. 


4.2  WEAK  SIGNAL  IN  NOISE  SUBOPTIMUM  PROCESSOR  REALIZATION 


Assuming  that 


S(w) j/N(u)  «  1  for  j  =1,  2,  ...,  J 
and  using  Equations  (3.5.1-17c-g),  Equation  (3.5.1-12)  becomes 


(4.2-la) 


|hj(«) 


SjM/NjU)  __  Sj(w)/Nj(u) 

1  +  Gj(u)  Sj(a>) /Nj(oj)  1  +  ”  SjW/W 


(4.2-lb) 


2  _  S.(oj)/N  .(co) 


a  .(a,)  =  Nj(u)  |hj(a>)  |  =  1  +“j  $  (<4/N(u>) 

\i 


(4.2-lc) 
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S^O))  Pi(co)  +  N(oj)  Q ( to) 

i=l 

ift 


s  Q-i(io)  Nj(M)/N(w) 


(4.2-ld) 


(4.2-le) 


u 

ij(«)  =  2  si(“)  +  N(u) 


(4.2-lf) 


Pj(w)  =  VjU)  1M  VJ(a>) 


S -(a)) /N  .(to)  r~  i  1 

bi<“)  *  1  ♦  th  trLQj  V"1  *i  VjH 

_  vl  _ 


(4.2-lg) 


(4.2-lh) 


and  the  simplified  likelihood  Equation  (3.5.1-12)  of  estimating  time  delay 
between  sensors  i  and  i  +  1  due  to  target  j  becomes: 


3A(t) 

W  * 


|  ju)  |h j (<*»)  a*(to)  Q'*(<o)  V j (oj)  [4>i  -  b'j(co)  1M] 


Vj(co)  Q"1  (to)  a(w)  -  T  b-j(u)  j 


(4.2-2) 


Note  that  Equation  (4.2-2)  results  in  a  substantial  simplification  of 
Equation  ('.5.1-12).  For  the  two-target  and  two-sensor  case  studied  in 
Section  3.,  ,  t.hf  esulting  two  likelihood  equations  (3.5.1  -17a)  and 

(3.5.1-17b)  become: 
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»A(tj  t2)  f  | 

- Fj -  =  J  j“  j  lhl !  a*  Q  v1($1  -  bj 


1„)  Vf  Q-1  o.  -  T  bjj  g 
(4.2 -3a) 


dco 


=  o 


3A(tl  t2) 


/  I 

V0 


=  0 


jo.  j  |h2|2  a*  Q"1  V2(^  -  b2  1M)  Q'1  a  -  T  b2j  ^ 

(4.2 -3b) 


where 


2  y(S2  +  N)‘ 


/N 


(4.2-3c) 


|h 


2  _  S2/(S1  +  *Y 


2 1  1  +  2  S2/N 


Q_i  =  I 


(4.2-3d) 

(4.2-3e) 


*1  = 


0  -1 

1  0 


(4.2-3f ) 


bl  * 


b2  = 


/  v  i 

/  VN  \  r 

\r+  2  s1  /N J 

U  +  2  S2/n;  [ 

(,  Y»  \ 

/  S2/N  \ 

\1  +  2  S^N  / 

\  i  +  2  s2/n; 

ju)(Tj-T2)  -jaiCrj-Tp) 

e  -  e  1  c 


(4. 2-3g ) 


JUiy-Ti)  -JOj(Tp-T,  ) 

e  -  e  c  1 


(4.2-3h) 


Using  Equations  (4.2-3c)  through  (4.2 -3h ) ,  Equations  (4.2-3a)  and  (4.2-3b) 
can  be  further  simplified  tc 
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9A(t1  t2) 


oo 

-hf 


SjSg/N" 


GXlx2  +  (r+  ?  Sj/n)  (i  +  2  s2/n) 

(l*j  |hj  vj  q'1  2!2)e'J“T2jeJ“Tl 


(4.2-4a) 


m(t1,  t2>  .  T 
3  T0 


=  ^  / 


Jto  ;  |h„|2  G 


21  “xix2 


sx  s2/ir 

T^WJT-r?  "S27n") 


=  o 


dw 


(4.2 -4b ) 


T  1  A 

where  x  =  (x,  x9)  =  Q“  a,  Gv  is  the  estimated  cross  power  spectrum,  and 

x  ^  x  ^  x  2 

Vp  y_2  are  the  time  delay  steering  vectors.  A  block  diagram  of  this  processor 
is  shown  in  Figure  4-1.  Note  that  given  two  sensors,  a  single  channel  GCC 
under  a  multiple  target  environment  is  known  to  be  biased.  The  coupling  shown 
in  Figure  4-1  provides  the  required  bias  correction.  Note  that  the  bias  cor¬ 
rection  term  is  a  function  of  SNR  and  INR.  When  the  power  spectra  are  not 
known,  the  quantities  SNR  and  INR  must  be  substituted  by  their  estimated 
values.  Thus,  the  optimum  power  spectral  estimator  discussed  in  Sec¬ 
tion  3.5.3  is  applicable.  Finally,  if  the  signal  spectrum  and  interference 
spectrum  are  separable  (i.e.,  no  overlapped  region),  the  two  channels  become 
uncoupled. 


4.3  A  SINGLE  TARGET  ASSUMPTION  SUBOPTIMUM  PROCESSOR 

For  completeness  we  now  include  a  study  on  the  suboptimum  multitarget 
processor  using  a  single  target  formulation.  The  resulting  processor  perform¬ 
ance  has  been  studied  widely  in  the  context  of  the  localization  variables 
(i.e.,  range  and  bearing)  (References  7,  10,  11,  12,  23,  24).  However,  a 
direct  study  on  the  time  delay  variable  has  not  been  seen  in  the  open 
literature. 
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A  single  target  processor  can  be  obtained  from  the  multi  tar  get  proces¬ 
sor  by  setting  the  interference  power  spectra  to  zero.  For  the  two-sensor 
case,  the  resulting  processor  reduces  to  the  GCC  processor  discussed  in 
Section  3.5.1.  The  time  delay  is  obtained  from  the  GCC  by  locating  the  peak 
of  the  GCC  function  (assuming  that  SNR  is  sufficiently  high  so  that  a  dominant 
peak  can  be  detected).  In  the  presence  of  interference  the  resulting  esti¬ 
mates  are  known  to  be  biased.  In  addition,  they  affect  the  time  delay  variance 
performance.  Here  we  shall  quantify  the  performance  in  more  detail. 

Consider  a  general  two- sensor  J  target  problem.  Let  the  true  time  delay 
to  target  j  be  t.  for  j  =  1,  2,  ...,  J.  Then  the  frequency  domain  representa- 
tion  from  T  seconds  of  observation  time  can  be  written  as 


alk  =  -lk  Bk  +  nlk 


“2k  =  -2k  Bk  +  n2k  ’  k  ■  1,  2,  ....  B 
where  the  complex  vectors  v1(<,  and  Bk  are  defined  by 


.J“k0ii 


*“k°12 


H  Du 


w) 


1,  2 


(4.3-la) 

(4.3-lb) 


(4.3-2a) 


and 


Sk  * 


*  ( 


Bkl  Bk2 


kJ 


)T 


(4. 3-2b) 


where  D.^  is  the  propagation  time  delay  from  target  j  to  sensor  i  and  Bkj-  is 
the  Fourier  component  of  the  signal  spectrum  of  target  j.  The  peak  of  the  GCC 
is  obtained  by  locating  the  null  of  the  function  (see  Equation  3.5.1-14c) 
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f(x) 


u 

•  £  n  in 


JWbT 


alk  a2k 


k=-B 


(4. 3-3a) 


2 

where  |hk|  ,  the  spectral  shaping  filter,  is  given  by 

lhkl  '  1  +  1  Skl/Nk  *  (4. 3 -3b ) 

Note  that  without  loss  of  generality  we  have  let  j  =  1  be  the  target  of  interest 
and  let  the  remaining  0-1  targets  be  interferences.  An  alternate  selection  of 
the  frequency  shaping  filter  is  the  multi  tar  get  spectral  shaping  filter  (see 
Equation  (3.5 . l-4b ) ) : 


Skl^kl 


1  +  Gkl  Skl/Nkl 


(4.3-4a) 


where  from  Equations  (3.5.1-2c)  and  (3. 5. 1-3) 


+  N 


and 


3kl  =  ^kl  Gkl  ^kl 


»1 


N 


kl (S  Skj  Pkj  +  Nk  Qk  )  1 
yj=2  / 


\\ 


(4. 3-4b) 
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~  -i 

'  v  V  hi 


(4.3-40 


for  the  two  sensors  with  =  I. 


Therefore,  a  simpler  form  of  Equation  (4. 3 -4a)  is 


Skl^kl 

nr\T7^ 


(4.3-4d) 


For  example,  for  the  two-target  case.  Equation  (4.3-4d)  reduces  to 


,2  .  Skl/(Sk2  +  V  . 

k  '  1  +  2  A 


(4. 3-4e) 


The  basic  derivation  of  the  bias  and  variance  is  shown  in  Appendix  H 
(Equations  (H-8e)  and  (H-17b)).  The  expressions  for  the  bias  and  the  variance 
for  the  two-target  case  with  identically  flat  signal,  interference,  and  noise 
power  spectra  are  given  by  Equations  (H-12)  and  (H-21a)  as 
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m 


and 


VAR(Xl)  =  — 
1  T 


1  -  (a^  p(2A^)  +  a 2  p(2A2)  +  2a^  p(A^  +  A2)) 


2( ax  p(A:)  +  a2  p(A2))2  R( 0) 


(4. 3 -5b ) 


where  S/I  is  the  si gnal-to-interf erence  ratio  and 


a,  = 


S/N 

1  +  S/N  +  I/N 


(4. 3-6a) 


a„  = 


I/N 


2  "  l'TT/N  +  I/N 


(4.3-6b) 


Ai  =  t,  -  T 


1  ll  ‘  ll 


(4.3-6c) 


A2  =  T1  *  t2 


(4. 3-6d) 


CO- 


R(t)  -l  cj  cos  wt  duj 


(4.3-6e) 


p(x)  =  R(t)/R(0)  . 


(4. 3-6f ) 


Note  that  in, the  limit  as  Aj  -►  0,  A?  -*•  ®,  we  have  p(A2)  -*■  0,  p(Aj)  -►  1  and 
the  steady  state  variance  is 
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VAR  JV 


Tim  VAR(tx) 
Ax  ♦  0 


Therefore,  one  can  define  a  normalized  variance  as 


(4.3-7) 


var(Ti) 


p(A^)  +  2a^  p(A^  +  A2)  +  p( Ag) ) 

2(ax  p(Aj)  +  a2  p(A2))2 


(4.3-8) 


Note  that  Equation  (4.3-7)  does  not  reduce  to  the  single  target  case 
because  of  the  presence  of  spectral  interference.  This  is  certainly  true  for 
the  case  of  two  omni-directional  sensor  arrays  since  interference  power  is  not 
spatially  attenuated  as  a  function  of  time  delay  separation.  However,  for 
sensor  arrays  with  large  base  length  separation,  this  condition  is  again  satis¬ 
fied  since  a  small  spatial  separation  could  produce  a  large  time  delay  separa¬ 
tion  and  insignificant  spatial  attenuation. 

4.4  GCC  PERFORMANCE  IN  THE  PRESENCE  OF  INTERFERENCE 

Performance  of  the  GCC  in  the  presence  of  interference  is  compared  to  the 
optimum  processor  in  the  remainder  of  this  section.  Figure  4-2  shows  the 
interference  of  the  expected  value  of  tie  GCC  due  to  signal  only  by  the  GCC 
from  a  second  target  with  an  identical  spectrum  shape  but  a  3  dB  smaller  signal 
power.  It  can  be  seen  that  (1)  the  two  GCCs  are  merged  to  one  (i.e.,  it  fails 
to  resolve  the  target  from  the  interference),  and  (2)  the  peak  of  the  combined 
GCC  is  biased.  Figure  4-3  shows  the  same  GCC  in  the  form  of  a  3-D  interfer¬ 
ence  pattern  as  a  function  of  target-interference  time  delay  separation. 
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Figure  4-2.  Multitarget  Time  Delay  Interference  (SNR  =  10  dB,  INR  =  7  dB) 


Multitarget  Time  Delay  Interference  as  a  Function  of  Time  Delay  Separation 
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(Note  the  combined  GCC  in  Figure  4-2  corresponds  to  the  curve  with  unity  sepa¬ 
ration  in  Figure  4-3.)  Since  the  MLE  is  asymptotically  an  unbiased  estimator, 
the  optimum  processor  has  no  bias  for  sufficient  integration  time.  Therefore, 
it  resolves  the  multitarget  ambiguity.  This  is  the  primary  advantage  of  the 
optimum  multitarget  processor.  Figure  4-4  shows  the  bias  characteristics  of 
the  GCC.  Note  that  when  the  time  delay  separation  is  small,  the  bias  is 
proportional  to  the  separation,  indicated  by  Equation  (4.3-5a).  Figure  4-5 
presents  the  resulting  GCC  variance.  The  GCC  variance  about  the  estimated 
mean  is  normalized  by  the  corresponding  variance  with  no  interference  (Equa¬ 
tion  (3.5. 1 -16b ) ) .  Also  shown  in  Figure  4-5  is  the  CRLB  of  the  optimum  two- 
target  processor.  This  optimum  processor  not  only  resolves  the  bias  but  also 
has  a  smaller  steady  state  variance.  Note  that  even  with  no  time  delay  inter¬ 
ference,  the  variance  of  the  GCC  does  not  approach  the  bound  because  of  the 
presence  of  interference  in  the  frequency  domain.  Also  note  that  with  no  time 
delay  separation  between  target  and  interference,  the  CRLB  is  singular, 
reflecting  the  inherent  inappropriateness  of  a  two-target  formulation  to  a 
one-target  estimation  problem. 

Also  shown  in  Figure  4-5  is  the  total  Root  Mean  Square  (RMS)  error  of  the 
conventional  GCC  processor  in  the  presence  of  interference.  Note  that  the 
CRLB  of  the  optimum  two -target  processor  results  in  a  significantly  reduced 
total  RMS  error.  Finally,  we  remark  that  identical  signal  and  interference 
power  spectra  are  the  worst  conditions.  For  different  signal  and  interference 
spectra,  the  optimum  processor  is  well  behaved  even  at  zero  time  delay  separa¬ 
tion  (see  Figure  3-7). 
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CHAPTER  5 

IMPROVED  MULTITARGET  PARAMETER  RESOLUTION 


5.1  INTRODUCTION 

In  Chapter  3  we  studied  the  optimum  multitarget  multisensor  parameter 
estimator  using  the  MLE  procedure.  In  Chapter  4  we  studied  a  number  of  sub¬ 
optimum  realizations.  We  pointed  out  that  the  GCC  processor  can  be  considered 
as  a  suboptimum  processor  in  a  multitarget  environment.  The  use  of  GCC  in  a 
multitarget  environment  results  in  a  poor  multitarget  parameter  resolution; 
i.e.,  it  yields  a  single  biased  estimate  when  separation  between  targets  is 
small.  The  use  of  an  Optimum  Multitarget  Processor  (OMP)  can  provide  a  sig¬ 
nificant  improvement  in  resolution  (as  shown  in  Figure  4-5).  Unfortunately, 
the  OMP  also  requires  a  major  modification  of  many  existing  systems  where  GCC 
processors  have  already  been  implemented.  Therefore  in  this  chapter  we  inves¬ 
tigate  alternate  procedures  for  improved  multitarget  parameter  resolution. 

The  particular  approach  to  be  studied  is  the  post  GCC  multitarget  pro¬ 
cessor.  In  this  approach,  additional  multitarget  processing  capability  is 
provided  at  the  GCC  outputs.  Thus  with  this  approach,  the  GCC  processor  needs 
no  modification.  In  effect,  the  existing  GCC  processor  can  serve  conveniently 
as  a  pre-processor  of  a  multitarget  estimator. 

This  chapter  is  organized  as  follows.  Section  5.2  derives  the  post  GCC 
multitarget  processor.  Section  5.3  presents  the  estimator  performance  bound. 
Section  5.4  discusses  the  simulation  procedure  and  results.  Specifically,  we 
compare  the  performance  of  the  conventional  GCC  processor,  the  post  GCC  multi - 
target  processor,  and  the  optimum  multitarget  processor. 


5.2  POST  GCC  MULTITARGET  PROCESSOR 

From  Equation  (3.5.1-14b),  the  GCC  output  can  be  written  as 
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R(t)  =  £ 
k=-B 


hk  I  alk  a2k  e 


H' 


(5.2-1) 


where  | | ,  and  are  defined  as  shown  in  Appendix  H.  Now  Equation 
(5.2-1)  can  be  written  as 


R(t)  =  R(t)  +  W(t)  (5.2-2a) 

where  R(t)  denotes  the  deterministic  component  and  W(x)  the  random  component 
of  the  noisy  GCC  ouput.  For  practical  application.  Equation  (5.2-1)  is  usually 
realized  via  a  Fast  Fourier  Transform  (FFT).  Thus  the  GCC  output  consists  of 
a  discrete  set  of  observations.  Therefore,  let  At  be  the  sampling  time,  then 
the  discrete  GCC  output  can  be  written  as 

R(nAt)  =  R(nAt)  +  W(nAt)  ;  n  »  0,  ±1,  ±2,  ....  ±N/2  .  (5.2-2b) 


Using  Equation  (5.2-1),  the  deterministic  and  the  random  components  of 
R(nAt)  can  be  obtained.  For  example,  taking  the  expected  value  on  both  sides 
of  Equation  (5.2-1),  the  deterministic  component  is 


B 

R(nAt)  =  ^ 
k=-B 


l\l  °lk  afk  e 


jw^nAt 


l5.2-3a) 


But  from  Appendix  I,  we  have  the  relation 
J 

Sk,  e  ^  (5.2-3b) 

1*1 

where  t.  is  the  true  time  delay  between  two  sensors  for  target  i  and  Ski  is  the 
corresponding  target  discrete  power  spectral  density.  Substituting  Equation 
( 5. 2 -3b)  into  Equation  (5.2-3a)  yields 
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G22(k)  =  =  S  Ski  + 


(5.2-4c) 


Gl2(k)  =  aik  =  X]  Ski  e  ^  1  • 

i=l 


(5.2-4d) 


Now  define  the  2J  unknown  parameter  vector  by 


9.  ( Sj ,  S^j  G-jJ  »  T^,  •••>  T-j) 


j,  ip  i2, 


(5.2-5a) 


and  the  assumed  matching  function  by 


V®)  3  S  Si  pi(nAt  ”  Ti)  » 

i=l 


( 5 . 2 -5b) 


then  the  observation  Equation  ( 5 .2-2b)  can  be  written  as 

R(nAt)  =  hn(e)  +  W(nAt)  ;  n  =  0,  ±1,  ±2,  . . . ,  ±N/2  .  (5.2-5c) 

In  matrix  notation,  this  can  be  written  as 


Z  =  h(e)  +  W 


(5.2-6a) 


where 


Z  =  [R{ -NAt/2) ,  R( NAt/2) ] 


(5 .  '.'-6b) 


jl(±)  ^_|\jy2(9. )»  •  ••»  ^  — N/ 2  ( — ^ 


(5.2-6c) 
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W  =  [W( -NAt/2),  ....  W( NAt/2) ] T  (5.2-6d) 

are  N+l  dimensional  vectors. 

Note  that  the  noise  vector  is  zero  mean  with  matrix  covariance 

E(W  WT)  =  A  (5.2-6e) 

where  the  mn  element  of  A  is  given  by  Equation  (5.2-4a).  Using  an  IMS  error 
criteria,  the  best  estimate  of  the  unknown  parameter  vector  is  obtained  by 
minimizing  the  function 

J(e)  =  [Z  -  h(i)lT  1 1  *  Mi)]  •  (5.2-7) 

Thus  the  best  estimate  of  e  is  given  by 


e  =  Min  Arg  J(e) 

9 

or  equivalently  the  vector  null  equation 


(5.2-8a) 


(5.2- 8b ) 


Note  that  because  we  are  trying  to  find  a  best  match  of  the  observed  GCC  out¬ 
put  to  an  assumed  reference  function,  the  resulting  estimator  is  appropri¬ 
ately  called  the  Matched  Estimator. 


5.3  ESTIMATOR  PERFORMANCE  EVALUATION 

In  this  section  we  present  the  multi -parameter  covariance  matrix  bound 
for  the  post  GCC  multitarget  processor.  Rewriting  the  vector  null  equation 
(5.2-8b)  as 
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y(e)  = 


3h(e) 

le- 


T 


2  -  h(e)j  =  0  , 


(5.3-1) 


the  Taylor  series  expansion  of  y(e)  about  e^,  the  true  parameter  vector  yields 
(ignoring  the  higher  order  terms) 


3 l(o) 

*<t>  *  £(»0) +  sr- 


(i  "  ®q)  • 
£  = 


(5.3-2a) 


Since  by  definition  y(e)=0,  we  must  have 


i<80> 


a l(») 


3e 


<t  - 
t  *  to 


=  -  A(©q)  Se 


(5.3-2b) 


^  A  ^ 

where  6e  -  (e  -  and 


^(®) 

A(e0)  =  3e 


£  =  £o 


Shje) 

3e 

[Z  -  h(e0)]  - 

' 9h (e ) ' 

3e 

T 

3h(e) 

3e 

£  = 

l<*> 

II 

£ 

£  =  £o 


~  ”H  ( ®q )  H(©q) 


(5.3-2C) 


In  deriving  Equation  (5.3-2c),  it  was  assumed  that  the  first  term  is  negligible 
for  sufficiently  high  SNR.  Therefore,  from  Equation  (5.3-2b),  one  obtains  the 
matrix  covariance  equation 
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COV [^(®q) ]  =  A(9q)  COV[6e]  A^)  .  (5.3-2d) 

But  from  Equation  (5.3-1) 

COV^Sq)]  =  HT(eQ)  A  HCsq)  .  (5.3-2e) 

Therefore,  the  estimator's  matrix  covariance  is 

COV [6e]  =  A’1(e0)  HT{©0)  A  H(e<))  A'1^)  .  (5.3-3) 


5.4  SIMULATION 

A  computer  program  was  developed  on  the  VAX-11/780  at  NUSC  to  simulate 
the  multitarget  GCC  output.  For  simplicity,  a  two-sensor,  two-target  envi¬ 
ronment  was  chosen.  The  GCC  output  observation  window  was  set  at  ten  times 
the  reciprocal  signal  bandwidth.  Both  targets  are  assumed  to  have  a  broadband 
flat  spectra.  One  thousand  Monte  Carlo  iterations  were  used  for  each  time 
delay  separation.  Simulation  results  are  plotted  along  with  the  theoretical 
performance  predictions. 


5.4.1  Simulation  Procedure 


Figure  5-1  diagrams  the  simulation  procedure.  The  multitarget  GCC 
output  observation  vector  (Equation  (5.2-6a))  was  generated  by  adding  observa¬ 
tion  noise  sequence  with  prescribed  covariance  to  the  noise-free  GCC  com¬ 
ponent.  The  noise  vector  was  generated  using  the  following  procedure: 

(1)  calculating  the  multitarget  GCC  covariance  matrix  from  Equation  (5.2-4a), 

(2)  factoring  the  matrix  into  lower  and  upper  triangular  matrices  using  the 
Gramm-Schmidt  orthogonalization  procedure,  and  (3)  multiplying  the  lower  trian¬ 
gular  matrix  by  a  white  noise  vector  to  produce  the  desire  correlated  observa¬ 
tion  noise  vector.  Finally,  the  GCC  observation  vector  was  processed  by  the 
matched  estimator. 
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5.4.2  Discussion  of  Results 


Since  a  detailed  investigation  of  the  matched  estimator  will  be  shown 
elsewhere  (References  27,  28),  only  a  brief  discussion  of  the  simulation 
results  will  be  presented  here.  Figure  5-2  shows  the  bias  characteristics  of 
a  conventional  6CC  estimator;  the  target  assumed  a  broadband,  one-sided 
spectrum  between  0  and  B  Hz  while  the  interference  assumed  a  broadband,  one- 
side  spectrum  between  B/2  and  B  Hz,  the  target  strength  is  at  0  dB  and  the 
interference  is  at  -1  dB.  Note  the  large  bias  when  targets  are  separated  just 
below  the  reciprocal  signal  bandwidth.  In  this  region,  targets  are  not 
resolved  with  the  conventional  GCC  estimator.  Figures  5-3  and  5-4  show  the 
normalized  rms  (normalized  by  the  time  delay  standard  deviation  of  target 
only)  performance  of  a  conventional  GCC  estimator,  the  post  GCC  matched  esti¬ 
mator,  and  the  optimum  multitarget  estimator  as  a  function  of  time  delay 
separation.  While  Figure  5-3  assumed  identical  broadband  signal  and  inter¬ 
ference  spectra.  Figure  5-4  assumed  the  interference  occupied  the  upper  half 
frequency  band  but  maintained  the  same  interference  power.  Note  that  the 
conventional  GCC  rms  was  clipped  at  a  degradation  ratio  of  20  for  time  delay 
separation  less  than  the  reciprocal  signal  bandwidth.  Both  Figure  5-3 
and  5-4  show  the  marked  improvement  of  the  matched  estimator  over  the  GCC 
estimator.  Note  the  close  performance  between  the  matched  estimator  and  the 
optimum  processor. 
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Figure  5-1.  Post  GCC  Multitarget  Processor  Simulation 
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Figure  5-2.  GCC  Bias  Versus  Time-Delay  Separation  (SNR  =  0  dB,  INR  = 
S : 0  +  8,  I: B/2  -  B,  N:0  -  B) 
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CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 


We  have  studied  in  detail  the  methodologies  of  optimum  signal  processing 
for  passive  time  delay  estimation  in  a  multisensor,  multitarget  environment. 
Our  investigations  were  motivated  by  (1)  the  inherent  inability  of  the  exist¬ 
ing  time  delay  estimator  to  resolve  estimation  bias  in  a  multitarget  environ¬ 
ment,  and  (2)  the  apparent- lack  of  research  and  understanding  in  this  area. 

In  the  literature,  the  studies  of  interference  were  confined  to  studying  the 
effect  of  interference  on  the  existing  processor  and  the  methods  of  inter¬ 
ference  suppression.  The  location  or  direction  of  interference  was  usually 
assumed  known.  Our  study  as  presented  here,  however,  treats  the  interference 
as  another  target  of  interest. 

We  pointed  out  that  the  traditional  approach  is  optimum  for  a  single 
target  only.  In  the  multitarget  environment,  the  existing  approach  is  biased 
due  to  the  inherent  mismatch  between  the  single  target  signal  processor  design 
and  the  multitarget  operating  environment.  We  argue  that  for  a  high  perform¬ 
ance  sonar  system,  it  is  important  to  have  an  unbiased  processor  since 
measurement  bias  cannot  be  easily  removed  with  further  post-estimator 
processing. 

The  optimum  multisensor,  multitarget  time  delay  processor  is  derived 
from  a  Maximum  Likelihood  viewpoint.  The  actual  processor  is  obtained  by 
reducing  the  likelihood  equation  via  straightforward  but  somewhat  tedious 
manipulations  to  the  simplest  form.  The  general  multisensor,  multitarget 
processor  is  applied  to  a  two-sensor,  one-target  case.  The  resulting  proces¬ 
sor  is  shown  to  be  identical  to  the  GCC  studied  by  Knapp  and  Carter  (Ref¬ 
erence  1).  We  next  studied  the  two-sensor,  two-target  problem.  We  showed  the 
optimum  two-sensor,  two-target  processor  and  presented  the  CRLB.  We  then 
extended  our  study  to  address  the  estimation  of  localization  parameters.  Much 
of  our  attention  had  been  spent  on  the  discussion  of  a  basic  three-sensor 
ranging  array.  We  showed  that  for  passive  ranging  and  directional  finding,  an 
optimum  processor  is  the  focus  beamformer  which  yields  a  direct  estimate  of 
range  and  bearing.  We  argued  that  an  alternate  approach  is  to  measure  the 
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inter-sensor  time  delay  and  then  geometrically  map  the  time  delay  measure¬ 
ments  to  the  corresponding  ranges  and  bearings. 

The  focus  beamformer  approach  can  be  called  the  one-step  approach, 
whereas  the  alternate  time  delay  approach  can  be  called  the  two-step  approach. 
It  was  shown  in  Section  3.5.2  that  both  approaches  yield  identical  perform¬ 
ance  in  terms  of  the  CRLB.  However,  there  are  major  differences  between  the 
two  approaches:  (1)  the  focus  beamformer  requires  searching  over  a  range/ 
bearing  space  of  a  correlation  function  which  is  asymmetric  with  respect  to 
the  range  and  bearing  variables;  on  the  other  hand,  any  correlation  over  the 
time  delay  variables  is  always  symmetrical;  (2)  for  tracking  purposes,  the 
focus  beamformer  approach  requires  a  two-dimensional  error  detector  design 
whereas  the  time  delay  approach  can  be  implemented  using  two  one-dimensional 
error  detectors;  and  (3)  arguing  from  the  Law  of  Large  Numbers,  both 
approaches  yield  Gaussian  measurement  noise.  However,  for  the  time  delay 
approach,  the  resulting  range  and  bearing  estimates  could  be  biased  and  have 
non-Gaussian  statistics  if  a  direct  non-linear  geometric  mapping  from  time 
delay  measurements  is  used.  For  a  practical  implementation,  the  two-step  time 
delay  approach  is  preferred  because  the  symmetry  of  the  correlation  function 
over  the  time  delay  variables  yields  simple  and  efficient  tracking  logic. 

This  is  especially  advantageous  for  a  low  SNR  environment  where  the  increased 
smoothing  time  required  {to  make  the  estimator  efficient)  can  be  achieved  via 
a  simple  feedback  design.  Finally,  for  the  three-sensor  ranging  array,  the 
potential  bias  and  non-Gaussian  statistics  can  be  minimized  if  mapping  to 
intermediate  variables  is  used  instead  of  range  and  bearing;  for  example, 
mapping  to  cosine  bearing  and  inverse  range  since  they  are  linearly  related  to 
time  delay  measurements. 

We  also  investigated  the  general  optimum  inter-sensor  time  delay  vector 
estimator.  Our  study  showed  that  given  M  sensors,  the  M-l  time  delay  can  be 
obtained  with  M-l  correlators.  This  is  in  marked  contrast  with  Hahn's 
approach  (Reference  25)  where  a  total  M(M-l)/2  correlations  are  required.  In 
addition,  we  presented  a  simple  expression  of  the  CRLB  of  time  delay  vector 
estimation  under  a  single  target  assumption.  Finally,  in  our  study  of  the 
three-sensor  ranging  array,  we  pointed  out  the  improved  performance  of  the 
optimum  formulation  versus  the  conventional  approach. 
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We  stated  that  one  of  the  strongest  assumptions  we  made  in  studying  the 
optimum  time  delay  processor  is  the  assumption  of  known  target  power  spectrum. 
Therefore,  we  briefly  addressed  the  problem  of  power  spectral  estimation.  We 
studied  the  optimum  spectral  estimator  for  the  two-sensor,  one-target  case 
and  the  two-sensor,  two-target  case.  We  briefly  studied  the  problem  of  joint 
time  delay  and  spectral  estimation.  A  somewhat  surprising  result  is  that  time 
delay  estimates  and  spectral  estimates  are  uncorrelated.  This  implies  that 
the  joint  time  del  ay /spectral  estimation  does  not  degrade  the  resulting  esti¬ 
mator  performance  when  assuming  either  is  known.  Furthermore,  we  found  that 
while  the  time  delay  CRLB  decreases  as  the  inverse  of  observation  time,  the 
power  spectral  CRLB  decreases  as  the  inverse  square  of  observation  time. 

We  next  addressed  the  problem  of  practical  implementation  of  the  optimum 
multisensor,  multitarget  time  delay  processor.  We  noted  that  one  approach  is 
to  assume  a  weak  signal  in  noise  environment.  Here  the  intricate  structure  of 
the  optimum  multisensor,  multitarget  time  delay  processor  is  drastically 
reduced  to  a  manageable  form.  We  also  discussed  a  single  target  assumption 
approach  as  a  suboptimal  processor  in  a  multitarget  environment.  We  provided 
a  numerical  illustration  of  the  single  target  processor  behavior  in  the 
presence  of  interference.  Performance  of  the  single  target  processor  was  com¬ 
pared  to  the  optimum  multitarget  processor.  We  noted  that  in  the  presence  of 
interference,  the  single  target  processor,  in  general,  was  biased  and  had  a 
larger  variance  except  when  the  target  interference  time  delay  separation  is 
small;  i.e.,  less  than  a  correlation  pulse  width.  Within  this  region,  how¬ 
ever,  the  optimum  multitarget  processor  has  a  variance  which  grows  without 
bound  as  the  separation  decreases  for  the  case  of  identical  signal  and  inter¬ 
ference  spectrum.  This  reflected  the  inappropriateness  of  using  a  multi¬ 
target  formulation  in  a  single  (merged)  target  environment. 

In  a  multisensor,  multitarget  environment,  an  optimum  processor  remains 
optimum  so  long  as  the  actual  number  of  sensors  and  targets  matches  the  number 
assumed  in  the  optimum  processor  design.  A  mismatch  in  either  the  number  of 
targets  (addressed  in  this  study)  or  the  number  of  sensors  (caused  perhaps  by 
element  failure)  will  automatically  degrade  the  processor  performance. 
Therefore,  a  key  element  in  using  the  optimum  multitarget  processor  is  the 
correct  detection  of  the  number  of  targets. 
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Last  but  not  least,  the  optimum  processors  which  we  have  derived  and 
studied  are  based  on  a  SPLOT  assumption.  For  a  moving  target,  this  assumption 
is  difficult  to  satisfy.  Therefore,  the  optimum  multitarget  processor 
presented  in  this  study  must  be  further  refined  to  account  for  the  moving 
target  environment. 
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APPENDIX  A 

DEFINITION  OF  COMPLEX  GAUSSIAN  PROBABILITY  DENSITY  FUNCTION  (pdf) 


This  appendix  briefly  describes  the  meaning  of  a  complex  Gaussian  proba¬ 
bility  density  function  (pdf).  A  thorough  treatment  of  this  subject  can  be 
found  in  Goodman  (Reference  26). 

Let  X(t)  be  a  zero  mean  wi de-sense  stationary  white  Gaussian  random 

o 

process  with  correlation  function  E£X(t)  X(t)]  =  a  5(t  -  t),  then  its  Fourier 
coefficients  from  a  T-second  observation  are  given  by 


T  / 


-jWt,T 

X(t)e  dt 


=  L  -  j  ;  k  =  1,  2,  ...,  B 


(A-l) 


where  1^  and  Qk  are  known  as  the  in-phase  and  quadrature  phase  components.  It 
can  be  easily  verified  that  Ik  and  Qk  are  Gaussian  distributed  with  the  fol¬ 
lowing  statistics: 


E(I.  )  =  E(QJ  =  E(  I.  Q.  )  =  0 


(A-2a) 


E(!k)  =  E(Qfc)  =  2T  * 


(A-2b) 


Let  =  ( Ik  Qk)  ,  then  the  bivariate  real  Gaussian  pdf  of  is  given  by 


p<4>  "f2")'1  V 


(A-3a) 


where 


\-£(i  iVl  0  c2 

K  H<  FT  Lo  0^ 


( A-3b) 


A-l 
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Now  a  complex  Gaussian  uni -variate  pdf  of  Xk  is  defined  by 


P<Xk>  =  "'1  I'ltl'1  ck_1  V 


where 


ck  =  E(Xk  Xk}  =  °  /T  ‘ 

It  is  straightforward  to  verify  that  IRJ"1  =  2|ck|"1  and  X*  = 

4T  Rj^  Z^/2,  hence  one  can  write 


P(Xk)  =  P(4)  =  (tto2)'1  exp 


r  2  ,  a2 


i _ II  \ 

°2  I 


Therefore,  in  general  let  _X  be  a  complex  Gaussian  B-dimensional  vector  such 
that 

X  =  I  -  j  3  (A- 


E(i  IT)  «  E(a  sT) =  V/2 

e(I5T)  =  -e(£1t)  =  -W/2 


Now  define  a  2B-dimensi onal  real  vector  as 


Z  =  (f  5T)T 


Then  it  can  be  shown  that  P(XJ  is  a  complex  Gaussian  pdf  defined  by 
P(X)  =  P(Z)  =  tt”B  |c j "1  exp{-X*  c-1  X} 
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APPENDIX  B 

CALCULATION  OF  CRAMER-RAO  LOWER  BOUND  FOR  TWO- SENSOR,  TWO-TARGET  CASE 

This  appendix  calculates  the  Cramer-Rao  Lower  Bound  (CRLB)  for  the 
two-sensor,  two-target  case.  Using  Equations  (3.4-1)  and  (3.4-5),  the 
CRLB  of  the  time  delay  estimates  are: 


VAR(?i )  >  [J-1]..  ;  i  =  1,  2  . 


( B-l ) 


Now  the  symmetric  Fisher's  Information  matrix  J  can  be  written  as 


■Jll 

r 

CNJ 

iH 

r-D 

r— 4 

II 

f — i 

1 

CNJ 

CNJ 

1 

-°i2" 

J21 

J22 

( J11 J22  "  J12} 

CNJ 

H 

1 

_ -J 

Jn_ 

Therefore 


VAR(^)  > - -r~  ;  i  =  1,  2  (B-2) 

1  "  (1  -  MJ2)  ii 


where 


M 


J12 


12 


( Jn  hz 


I1* 


is  defined  as  the  coefficient  of  mutual  dependence.  The  quantities  J-  . 

^  J 

are  defined  by  Equation  (3.4-5)  as: 


B-l 


1^1  1^  fa 
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aR-1  3Rk 


) 


i  =  1,  2 
j  =  1,  2 


(B-3) 


In  what  follows,  we  present  a  detailed  calculation  of  the  quantity  J...  Using 

^  J 

the  relations  given  in  Equations  (3.2-6)  and  (3.5.1 -4a) ;  i.e.. 


kl 


+  S 


k2  k2 


+  N,.  I 


( B-4) 


and 


1  «  1,  2 


(B-5) 


we  obtain 


3t,  3Tj  ) 


Skj  lhkil 


3Pk1  -  ,  3 Pi 
9Ti 


i;1 

3t  .  ) 

J 


3G, 


-  a 


'ki 


ki 


3Ti 


tr 


(\1  pkA'i  If) 


(B-6) 


The  trace  of  the  quantities  inside  the  parentheses  can  be  evaluated  as 
follows.  From  Equation  (3.2-7c),  we  find 


5k‘i  ■  Skl<Sk2  Pk2  +  Nk  D'1 

_  \l  fT  Sk2/Nk 
'Nk  i1  '  1  +  2  5k2/Nk  Pk2 


B-2 
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1 


‘  ak2 


~ak2  e 


J“lcT2 


"ak2 


e'Ja)kT2 


1 


"  ak2 


where 


Sk2/Nk 


*k2  "  r+-2  $ 


k2 


( B-7 


(b-s: 


and  similarly 


^k2 


k2 


1  -  a 


•nxi 


kl 


kl 


*akl  e 


JVi 


1  -  a 


kl 


( B-9] 


where 


Wk 

',1=rTT¥\ 


( B-10) 


Note  that  is  a  Hermitian  matrix  with  identical  diagonal  elements. 
Therefore,  letting 


TR  6757 


we  have 


-  3P.  •  .  qll  q12 

Oki  w1  -  n  .* 

°  q12  qll 


gj2 

qll  qJVj 


■Htj 


,Ja)kTj 


■>u  ’v‘ 


( B-ll) 


and  after  some  algebraic  manipulation,  we  obtain  from  Equation  (B-ll) 


5-1  !!jll  6-1  3P|<j  - 

gki  9t.  Mki  ' 


(B-12) 


where 


,•  2  jw.  (t  .+t  .) 

(q12>  e  J 


i  2  -JX(Ti-Ti) 
(qij)  e  ^  1  J 


A’J  =  (qjj)  [,{2  e>k(T,*T-')  -,*2  e‘J“k(TrTJ);  . 


(  «  «  £f)  -  - 


4  (Ajj  ♦  Ajf ) 


■2«j;  Re  {Ajj} 


(B-13) 


where  Re  {  }  denotes  the  real  part. 


B-4 
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Now  using  Equations  (8-7)  and  (8-8),  we  obtain 


tr 


21 


/5-l3Pkl\  -,2/\l\2 

[rn^r) 


(1  -  ak2) 


1  - 


\2 


(i  -  ak2) 


2*  cos  2cok A^ 2 


(B-14a) 


tr 


tn-lUkV  -22/\2\2  n 

\Qk2  3t2  /  2u)k\Nk  j  (1  "  akl^ 


21 


Nk?\  ^ 
k 


Jk2 


(1  -  akl) 


Y  cos  2cokA^2 


(B-14b) 


tr/5-1  5’1  3Pk2l  ?2/Nkirn  , 

1  \Qkl  aTl  Qkl  TET"]-  2wk  \  "NT  /  {1  ■  2a| 


k2)  cos  ukA12 


(B-14c) 


where 


t2 


is  the  time  delay  separation. 


On  the  other  hand,  we  have 


TR  6757 


i  JVi 
q12e 

i  _Ja)kTi 

q^e  k  1  +  q 

ni  pJVi 

qlle 

i*  -JATi 

qj2e  +  Q 

’ll 


’12 


i  * 


12 


’ll 


Thus  combining  Equations  (B-ll)  and  (B-15)  we  obtain 


tr 


3P, 


(St  pkt  Si  jp)  *  J 


n 


tr  < 


’ll 


’12 


1  * 


’12 


’ll 


q12 


qll 


( B11  q12  +  B12  qll}  e  k 


=  JVCij  -  cijl 

=  -2cok  Im{Ci  j  } 


where 


Cij  ^  (B11  q12  +  B12  qll)  e  ^  J 


•  2  .  ,_i  ,2  .noi 


■<qii>  e 


1J  j.  /„!  \  .  o  i 

+  (^12'  e  +  2  ^11  e 
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Im  {  }  denotes  the  imaginary  part. 


A-.At.-t-:  a- •  A  t.  +  t  . 

ij  -  1  j  *  ij  -  l  j 


Now  from  Equations  (8-7)  and  (B-8),  the  following  are  obtained: 


sr  f1  -  ak 2>2  1 


ak2  J"k412 
ak2 


(B-18a) 


N~ )  Cl  -  «ki>2  ( 1 


akl  JwkA21 

'T1kT 


(B-18b) 


12  "  V  N, 


(1  -  a„9r  U  - 


ak2  HA12\  "^12 

"  ak2  / 


(B-18c) 


Now  substituting  Equations  (B-18a-c)  into  (B-16)  yields 


trk\  pki  tt 


~-l  3Pkl\  _  ?  /\l  f  n  , 
qkl  3tj  J  ~  2uk\Nk  )  ak2  (1  ~  ak2^ 


2  sin  Vl2  “  T~ri^  sin  HA12 


( B-l 9a ) 


B-7 
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tr 


(1 


-  akl> 


2  sin  cO|^  A^2 


ak2 


k2 


sin 


and 


tr 


(5k‘  i 


~-l  aPk2 
pkl  «kl  3 TT 


2 a^ ^ ^  sin  * 


In  addition,  from  Equations  (3.5.1 -5b)  and  (B-ll),  we  have 


HTi 


i*  ■^a)kTi) 
q12  e  / 


and  for  i  =  1 


^kl  /  \l\ 

arj-  =  2“k  ak2  \  N^~ /  s1n  “kA12  * 


Finally,  substituting  Equations  (B-14a-c),  (B-18a-c)  and  (B-19a-c)  into 
(B-6)  we  find: 
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tr 


(.5L3a\- 

\  3x^  / 


Skl  lhkl  I  tr 


(<*£) 


-l  tr  /  J-l  p  O’1  3Pkl 

3kl  3tj_  1  \Qkl  Pkl  Qkl  3^ 


2\  Skl 


i-W2(H 


»  -  \2'  n 


where 


Yl  =  l 


=  jl  -Ei 


2 

cos  2(0^2  -  E2  sinc 


+  E3  (sin  a)kA12)  (sin  20)^)  J 


with 


Ei  = 


k2 


(1  -  ak2) 


E2  =  4  akl  I1  -  ak2> 


ft1)', 


E,  = 


2  akl  ak2 


(W 


(1  •  ak2'  k 


)| 

( B-22a) 

( B-22b) 

( B-22c) 

(B-22d) 


( B-22e) 


i* 


TR  6757 


m 


Similarly,  we  have 


(-#SKS2  u-*u)S 


tr 


L  2  U  U  2 


where  y2  is  obtained  by  exchanging  indices  between  1  and  2  from  y^. 
the  cross  term  is 


tr 


I 

\"  3-r,  ax2  / 


Sk2  |nkil2{f 


sPh  ~  i  aP 


k2 


3t-- 


/"-l  ~-l  9^k2  \ 

'  aki  3rf  tr(Qki  pki  Qki  a 


2{X  Sk2  lhkl  1 2  (1  “  2  ak2) 


(M 


Y12 


where 


Y12 


cos  u)j^  2  ~  2  i  a*,  o  l  hi  /  sin  uv  A 


kl  k2 


.2 


k  12 


Thus,  the  elements  of  the  Fisher's  Information  matrix  are: 


J11  =  2  S  \  Skl  lhkl  (n“)  (1  "  ak2)2  Y1 


k=l 


8  /N  \2 

J12  =  2  ^wk  Sk2  1  hkl  1 2  \N“)  (1  ‘  2ak2^  Y12 

k=l  K 


( B-23) 

Final ly, 

I 

I 

(B-24a) 

( B-24b) 

(B-25a) 

( B-25b) 


B-10 
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U  -  akl)2  Y2  • 


Using  the  definition  of  | hk i  p ,  Nki  and  aki  in  Equations  (3.5.1 -4b) ; 
(B-8),  and  (B-10),  Equations  (B-25a-c)  can  be  expressed  in  integral 
for  a  sufficiently  long  observation  time  T  as  follows: 


J11 


t  Fit  \( 1  *  vn  \ 

"  Tr  J  w  \1  +  G1S1/(  S2  +  N)  ) \l  +  2  S2/N /  Y1  da) 


12 


.  T  /  2  /  W"*  \/  1  \ 

*  IT  J  u  yl  +  61S1/(  S2  +_nT  )  \Y'+~? . 'S2'/N )  Y12 


da: 


and 


S?/N2 


22 


T  C  2  l  b2/IN"  \l  1  +  Sl/N  . 

"  TT Jr  “  ^1  +  G2  S2/(S1  +  N)  Jl  1  +  2  Sx/N  I  Y2  da) 


where 


Gx  =  (1  +  S2/N) 


/  S?/N  \  2 

2 -(n-hpr)  l1  +  6  12l 


G2  •  (1  +  S1/N) 


2  (  Sl/N  \  1 1  |  rJu)A12 ,2 

2  "^1  +  2  ^/N  J  I1  e  1 


( B-25c) 

(3.5.1 -2c ) , 
form 

( B-26a) 

( B-26b) 

(B-26c) 

( B-26d) 

( B-26e) 


B-ll 
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A  somewhat  simpler  expression  for  y^,  y^2,  and  Y2  can  be  obtained  as  follows. 
From  Equation  (B-22a-e)  we  obtain 


E 


1 


(1 


k2; 


(  VN  \2 

\1  +  2  S27Ny 

(  1  +  S2/N  ^2 

\r+  2  s27nJ 


s2/n  ^2 

l  +  S2/N/ 


(B-27a) 


E2  =  4  akl 


(i  -  ak2) 


( B-27b) 


where 


Sx/N  1 

'  1  +  s2/n 

1  +  G]_  S[/[S2  +  N) 

1  +  2  s2/n 

(B-27c) 


and 


E 


3 


2  akl  ak2  l\l\ 

(TTT^W  ) 


=  W2  Ex 


(B-27d) 


where 


Sl/N  -| 

s2/n 

[T  +  Gx  5^(52  +  N) 

1  +  2  S?/N 
/  ^ 

( B-27e) 
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Hence  Equation  (B-22b)  can  be  written  as 


=  1  -  ^cos  2(d)<A12  +  sin2  oj^A^  -  W2  sin  oj^A^  sin  2u>kA^2  j 
3 

=  1'ElEAn  Cosn  Vl2  (B"28a) 


where 


A0  =  4  1  +  Gx  sJ/(S2  +  N)  1  +T  S2/N 


1  +  s2/n 


Sl/N 

Ai  =  "4  i  +  Gj_  +  n; 


+  2  S2/N 


( B-28b) 


A2  =  1  -  aq 


A3  =  ~A1  * 


Similarly,  one  obtains 


3 

n=0 


^2  =  1  -  E2  Lj  Bn  C0S  WkA12 


(B-28c) 


where 


B-13 
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Finally,  from  Equation  (B-24b) 


Y12  =  cos  ^^12  ”  ^ 
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APPENDIX  C 

CRAMER-RAO  LOWER  BOUND  OF  TIME  DELAY  ESTIMATION  FROM  MULTI  SENSOR  ARRAY 


In  general,  the  Cramer -Rao  Lower  Bound  (CRLB)  is  given  by 


VAR  (e.)  > 


(C-l) 


where  J  is  the  Fisher  Information  matrix  whose  ij  element  is  defined  by 
(Equation  (3.4-5)). 


Under  the  assumptions  of  (1)  single  target,  (2)  spatially  incoherent  noise 
processes  and  (3)  identical  sensor  array  noise  power  spectrum,  the  CRLB 
can  be  evaluated  easily. 


From  Equations  (3. 5. 1-2)  and  (3. 5. 1-4)  we  obtain 


R 


k 


P,  +  N, 


( C-3) 


where  Nk  Qk  =  NkI  using  assumptions  (2)  and  (3),  and  the  relation 


C-l 
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sr:1  ,  3Pk 

k  •  M2  k 


3e.  '"k1  3ei 


( C-4) 


where 


ihi 


Sk/Nk 

1  +  M  VNk 


(C-5) 


and  M  is  the  number  of  available  sensors. 


From  Equations  (3.5.1 -10b )  one  obtains 


157 =  n  \  Ur  \ 


(C-6) 


where  V^,  the  steering  matrix,  is  given  by 


Vk.d1ag}eJ^eJ^...e^°M| 


( C-7) 


and  the  mn  element  of  the  matrix  J£-  is 

30  i 


(«r) 

mn 


_  3 


■g-  ( D  -  D  ) 
30^  m  n 


( C-8) 


Now  since 


aRk  .  s  9Pk 

30,  ‘  k  3e 
J  J 


(C-9) 


from  Equation  (C-3)  and  using  Equation  (C-4)  in  Equation  (C-2),  one  obtains 


C-2 
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the  expression 


[ 

\3«i  s»j  r 


\  ihki  tr 


\  3ei  3ej  / 


•sk  ihk  t2  ^ tr  ( vk  H-  <  vk  IS-  < )  •  <c-io> 


Now  using  the  relation  Vk  =  I  and  tr(AB)  =  tr(BA),  Equation  (C-10)  becomes 


'  J  *  J 


(C-ll) 


Therefore 


3k  l"k 


|  4  tr 


(39  39  \ 

ileTWrj 


B  c2/n2 

/  39  39  .2  Sk/Nk 

l^e-^j/Z-r^k  IT" M  Sk/Nk 


(C-12) 


Knowing  J.^'s  for  all  i  and  j.  Equation  (C-l)  can  be  used  to  calculate 

the  CRIB.  Using  Equations  (C-l)  and  (C-12),  we  shall  establish  the  following 

re lation. 


Given  M  sensors  under  a  single  target,  spatially  incoherent  and  spectrally 
identically  distributed  noise  environment,  the  CRLB  for  the  incremental  time 
delay  are  identical  and  equal  to 


C-3 
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VAR(xi)  > 


B 

k=l 


<4 


Vwk 


k  \1  +  M  Sk/Nk/ 


-1 


;  i  =  1,  2,  ....  M-l  .  (C-13a) 


Furthermore,  denote  the  time  delay  between  any  two  sensors  by  D,  then 


VAR(Dml)  »  VAR( ri )  for  all  i. 


(C-13b) 


Thus,  Equation  (C-13a)  implies  that  under  optimal  signal  processing, 
the  time  delay  estimation  between  any  two  sensors  improves  when  the  number 
of  sensors  is  increased. 


We  shall  establish  Equations  (C-13a)  and  (C-13b)  through  the  following 
intermediate  steps: 


(1)  Let  xi, 
vector,  then 


Tj  be  any 


two  elements  of  the  incremental  time  delay 


/a$_  3$_\  _  j"2  1  (M"J)  ;  J  -  1 
\ 3ti  dT.  J  j  _2  j  ( M-i )  ;  i  >  j 


( C-14) 


(2) 

VARU.)  > 


Sk/Nk 

1  +  M  VNk 


-1 


for  i  =  1,  2, 


M-l  .  (C-15) 


(3)  Define  D«  D  -  Dn,  and  x.  =  D.+, 


-  D . ;  then 
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VAR(Dmn)  =  VAR(t^)  for  all  m,  n,  and  i, 


( C-16) 


From  Equation  (C-8),  the  mn  element  of  the  matrix  $  is  given  by 


■  <Dmn>  ' 


( C— 17 ) 


Therefore,  the  matrix  $  can  be  written  as 


■(t1+t2) 


i  i  j. 

Ex* 

1-1 

M-2 


( t1+t2 )  t2 


0  .  . 


.  ( C-18) 


M-l  M-2 

Exi  Exi 

1-1  i=l 


Note  that  $  is  a  skew  symmetric  matrix  (i.e.,  <J>  =  -$  ) ;  all  elements  in  the 
main  diagonal  are  zeros;  the  incremental  time  delays  occupy  the  f.rst 
diagonal;  and  that  any  elements  above  the  first  diagonal  can  be  obtained  by 
summing  the  elements  to  the  left  and  below. 
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(C-20b) 


( C-20c) 


Thus,  using  Equations  (C-20b)  and  (C-20c)  in  Equation  (C-20a)  and  carrying 
out  the  necessary  matrix  operation,  one  finds 


C-6 
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tr  (H“  Ifr) =  -2i("-;i)  iJi' 

and  by  symmetry  argument,  we  obtain 

tr(l?75^)=  -2J(n,-i) 

This  proves  Equation  (C-14). 

Next  we  want  to  show  Equation  (C-15). 

Combining  Equations  (C-2),  (C-5),  (C-12)  and  (C-14),  we  obtain 


where  A  is  a  matrix  given  by 


( C-21a) 


( C-21b) 


(C-22) 


( M— 1 )  (M-2)  ...  1 
(M-2)  2( M-2)  ...  2 


TR  6757 


HI 


» 


r# 


■  i^'i  iw 


Thus  the  inverse  of  J  is 


J"1  =  A"1 


8  /  c2/m2 

2  /  Sk^Nk 

2  Z-f  “k  l  1  +  M  Sk/N. 
k=l  '  K  K 


-.-1 


(C-24a) 


-1 


and  therefore  the  diagonal  elements  of  J  are 


[J-1],,  ■  [A-1!,, 


r  B 

,  V'  2 

/  iK  \ 

2  2-4  wk 

k=l 

\TTirV\l 

i  *  1,  2,  . ..,  M-i  . 

( C-24b) 


By  inspection,  the  matrix  A-'*'  is  given  by  the  M  x  M  matrix 


( C-25a) 


or  the  ij  element  is  given  by 


2 

M 

1 


Ai/  =  V  M 


1  =  J 

I  i-J  I  -  1  • 

I  i  -  J  I  >  1 


( C-25b) 


-In  = 


This  can  be  verified  by  a  direct  multiplication  to  show  that  A-AA  =  I.  There¬ 
fore,  using  Equation  (C-25b)  in  Equations  (C-24b)  and  (C-15)  yields  the 
desired  result: 


C-8 
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VAR  (e.)  > 


B  /  .2  /n2  v: 

-  “k  (l  *  M  yNk  ) 


( C-26) 


This  proves  Equation  (C-15), 


It  is  easy  to  see  that 


COV  ( Ti , Tj )  = 


-VAR  (xi)/2  ;  if  | i - j |  =  1 

o  ;  if  I i - j I  >  1 


(C-27) 


Finally  writing 


mi  a 

■S’. 


( C-28) 


as  the  time  delay  between  sensors  n  and  m,  the  maximum  likelihood  estimation 

A 

(MLE)  of  D„„  is  linearly  related  to  the  MLE  of  t4.  Hence,  the  variance  of  D  „ 
mn  l  mn 

is  given  by 


m-1  m-1 

A  ^  a  A 

=  L  1  E[Ti  v  • 

i j  =n 


(C-29) 


C-9 


Using  Equations  (C-26)  and  (C-27)  in  (C-19)  yields 


VARfD^)  =  (m  -  n)  VAR(^)  -  (m  -  n  -  1)  VAR(^) 


=  VAR(Ti) 


Thus  the  optimum  MLE  of  time  delay  between  any  two  sensors  yields  the 
Cramer-Rao  Lower  Bound. 
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APPENDIX  D 

FUNCTIONAL  RELATIONSHIPS  BETWEEN  TARGET  LOCATION  VECTOR  AND 
MEASURED  TIME  DELAY  VECTOR 


This  appendix  develops  the  two-dimensional  mathematical  relationships 
between  target  location  parameters  and  time  delay  parameters  of  a  general 
three-sensor  array. 

Figure  D-l  shows  the  general  array/target  geometry.  In  principle, 
measurement  of  time  delays  between  sensors  Aj,  and  A2*  A3  Provldes  the 
necessary  set  of  relations  to  obtain  target  range  and  bearing.  Using  the 
Law  of  Cosine,  these  relationships  can  be  obtained  as 


Figure  D-l.  A  General  Three-Sensor  Passive  Ranging  Array  System 


D-l 
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T1  =  °2  ‘  D1 


I r  +  +  2rL^  cos  e 


(D-la) 


Jr^+L^  -  2rl_2  cos(e  -  <{>) 


( D-lb) 


where  D^,  i  =1,  2  and  3  are  the  propagation  time  delays  and  <J>  is  the  offset 
angle  of  one  sensor  with  respect  to  the  baseline  formed  by  the  remaining 
two  sensors. 

For  the  case  where  r  »  L^,  and  r  »  l_2,  Equations  (D-la)  and  (D-lb) 
become 


COS  9 


a  -  —  COS  9 


2  2 
sin^  9 


( D-2a) 


and  similarly 


2  .  2 

Lp  L2  sin  e 

X2  -  -  -  COS  9  +  — ^ - 


( D-2b) 


However,  the  exact  expressions  for  the  location  parameters  as  a  function 
of  time  delays  can  be  found  using  Equations  (D-la)  and  (D-lb)  as  follows. 


D-2 
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After  rearranging,  squaring  and  simplifying,  Equations  (D-la)  and 
(0-lb)  reduce  to 

|  (ctj)2  -  L^2  =  2r(ct1  +  Lj  cos  e)  (D-3a) 

|  (ct2)2  -  Lg1  =  -2r(cx2  +  l_2  cos(e  -  <f>)  .  (D-3b) 

We  first  solve  for  e,  the  bearing  and  then  solve  for  r,  the  range.  Dividing 
Equations  (D-3a)  and  (D-3b)  yields 


2  2 

ct^  +  cos  e  (ct^)  - 

ct2  +  L2  cos  (e  -  40  =  -  (CT2,2  .  l2  =  'k 

so  one  can  write 

Lj  cos  e  +  KL2  cos(e  -  <j>)  =  -(ct^  +  Kct2) 
and  from  which  one  obtains 

_1  /  CTi  +  Kct2  \ 

e  =  C  +  cos  (--4 - Z~) 


( D-4a) 


( D-4b) 


( D-5a) 


where 


A  sin  £  =  KL2  sin  4> 

A  cos  5  =  1-,+  KL,  cos  <b  .  (D-5b) 
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Therefore,  from  Equation  (D-5a),  one  obtains 


L  .  -1  /  ”1  +  Kct2  \ 

cos 

+  cos  \ - A - )\ 

(CT1 

+  Kct2)  (L1  +  KL2  cos  <{>) 

^A2  -  ( cx^  +  Kct2)2 

A2 

A2 

where 


A2  =  (L1  +  KL2  cos  <t>)2  +  (KL2  sin  4>)2 
=  (4  +  KL2)2  (1  -  a) 

where 

2K44  (1  -  cos  41) 

Ct  "  o  • 

(4  +  KL2)^ 

Now  substituting  Equation  (D-5c)  into  (D-3a),  one  obtains  the  range 
as 


r  = 


( C4 )  "  4 

2(ct-,  +  Li  cos  e) 


[(c^)2  -  L2]  (4  +  KL2)  (1  -  a) 

2{ct^(L^  +  KL2)(1  -  a)  -  +  Kct2)(1  -  6)  -  a  -  y2  KL 


D-4 


KL2  sin 
( D-5c) 

( D-5d) 

( D-5e) 

solution 


L2  sin  <{>} 


where 


8  = 


KL2(1  -  cos  <t>) 


Li  +  KC 


( D-6b) 


CT^  +  KCT2 

T  =  4*  KL2 


(D-6c) 


From  the  relation  (D-4a) 


■fc) 


2  1 


(?)' 


-(?) 


(D-6d) 


Equation  (D-6)  can  be  simplified  to 


r  = 


-(£>  -4  ■•(?)] 


(?)•(?)  •■<« 


(1  -  a) 


( D-7  a) 


where 


CTi 


r(<t>)  «  (a  -  8)  + 


&)..e),].vrw 


sin  <{> 


( D— 7b) 


is  a  function  of  the  offset  angle  <J>. 


Now  further  algebraic  manipulation  shows  that 
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Thus  Equation  (D-7b)  becomes 


CTi  -  KCT~  I  o 

r{<|>)  =  •  L  ~kl~  (1  -  cos  <*>)  +  “Vl  -  ot  -  y  sin  4> 


( D-7f ) 


-i  ‘'"a 


Note  for  a  co linear  array  system  (i.e.,  <j>  =  0),  we  have  a  =  0,  r(0)  =  0, 
and  the  range  and  bearing  equations  are  given  by 


Li  1  - 


•(?)’! 

*[(?)-(?)! 


( D-8a) 


■1  /  CT1  +  Kct2  \ 
V  L1  +  KL2  / 


e  =  cos 


(D-8b) 


In  addition  for  =  L2  =  L  and  r  »  L^,  using  Equations  (D-2a)  and  (D-2b),  it 
can  be  shown  that 


e  =  cos 


+  t2) 


sin^  e 

c<t2  -  Tl> 


(D-9a) 


(D-9b) 


If  e  is  defined  w.r.t.  the  broadside  direction.  Equations  (D-9a)  and  (D-9b) 
can  be  written  as 
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APPENDIX  E 

STATISTICAL  CORRELATION  BETWEEN  TIME  DELAY  ESTIMATES 
WITH  COMMON  INPUT  CHANNEL 


If  time  delay  estimates  are  obtained  from  two  Generalized  Cross- 
Correlators  (GCC)  whose  inputs  contain  a  common  noise  channel  (Figure  E-l), 
then  the  resulting  time  delay  estimates  are  correlated.  This  appendix  calcu¬ 
lates  the  resulting  correlation  from  a  frequency  domain  approach. 


Let  the  frequency  domain  (Fourier)  representation  of  a  T-second  observa¬ 
tion  of  waveforms  from  sensor  array  A^,  A2,  and  A^  with  signal  propagation 
delays  D^,  D2  and  D^,  respectively,  be  written  as: 


+  nlk 


(E-la) 


Ja)kD2 

“2k  '  6k  e  *  "2k 


(E-lb) 


a3k  =  Bk  e  +  n3k 


k  *  1,  2,  ...,  B 


(E-lc) 


with  Sk  =  E(BkB£)  and  Nik  =  E(nikn*k)  for  all  i  as  the  discrete  signal 
and  noise  power  spectra,  respectively,  at  frequency  =  2irk/T.  For  convenience 
we  assume  that  interarray  noise  processes  are  uncorrelated.  The  time  delay 
Tj  =  D2  -  D^  and  t2  =  02  -  ^  can  be  obtained  by  seeking  the  null  of  the 
GCC  functions  (Equation  (3.5.1 -14c ) ) . 


B 


fi(Ti>  ■  E 

k=-B 


★  * 

jwk  hlk  h2k  olk  a2k  e 


E-l 


Figure  E-l.  Generalized  Cross  Correlator  Pair  with  Common  Input  Channel 
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>(t2)  =  j 


H  h2k  h3k  ^k  “3k  e 


Ht2 


(E-2b) 


where  8  is  the  uppermost  frequency  of  either  the  signal  or  the  noise  and 
h^  is  the  frequency  shaping  filter  whose  frequency  response  is 


v4 


Vnik 


(E-2c) 


Let  t,  and  be  the  estimated  time  delay,  then  the  quantity  of  interest 
is  the  covariance  between  ^  and  t2;  i.e.,  C0V(t1,t2).  The  Taylor  series 
expansion  of  f ^ )  and  f2(r2)  about  the  true  time  delay  t°  and  t2  yields 


T?)  + 


( E-3a) 


f?(To)  - 


dfp ( To) 

W  +  -4f- 


xj)  +  .  . 


(E-3b) 


E-3 


m3* 
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Neglecting  the  higher  order  terms  (small  error  assunption)  and  setting 
f^l ( T1 )  =  ^2 ^ t2 ^  =  0  Equations  (E-3a)  and  (E-3b),  one  obtains 


( E-4a) 


3f2(^) 


( E-4b) 


We  shall  make  the  following  assumptions  as  a  result  of  a  long  observation 
time  process:  (1)  the  expected  value  of  the  derivative  is  equal  to  the  deriva¬ 
tive  of  the  expected  value,  and  (2)  the  derivative  is  uncorrelated  with  time 
delay.  Then  from  Equations  (E-4a)  and  (E-4b),  one  obtains 


3f 


St: 


(ti  - 


t?) 


(E-5a) 


(E-5b) 


where  for  simplicity  we  denote  the  expected  values  of  ^(1°)  and  f^x^) 
by  f°  and  f^.  Thus,  the  product  of  the  means  is 


—  —  3f?  3f° 


f°  -  1  2  ,  o  w  o> 

fl  f2  "  3^  (T1  '  Tl)  (t2  "  t2^  * 


(E-6) 


E-4 


* 


it 


Jtmm 


A 


A. 
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On  the  other  hand,  the  mean  of  the  product  is 


3fl3f!b - S77 - T7 

fl  f2  =  FJ  FJ  (T1  ‘  V  (t2  "  x2)  * 


(E-7) 


Therefore,  the  covariance  is 


C0V(t1,t2)  =  (tj  -  t°)(t2  -  t2)  -  (tj  -  x°)  (t2  -  T°) 


4:0  rO  ^0  ^0 

fl  '2  ■  '2 


3f  J  3f2 
3t1  3t2 


(E-8) 


From  Equations  (E-2a)  and  (E-2b),  the  quantities  in  Equation  (E-8)  can 
be  evaluated  as  follows. 


3fl 

fT  ’  Z<  (j“k> 


2  hlk  h2k  Sk 


(E-9a) 


3f2  ,  .  ,2 

FT  =  2^  ^“k^  h2k  h3k  Sk 


(E-9b) 


E-5 
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_  B 

fl  =  ^  Jwk  hlk  h2k  Sk 
k=-B 

_  B 

f2  =  ^2  J’“k  h2k  h3k  Sk 
k=-B 


and 

B  _ 

hi)<  h2k  h22,  h3$,  alk  “2k  “22. 

2,=-B 


f°  f°  - 
fl  r2 


k=-B 


But 


*  * 


“lk  “2k  “21  “31  ‘  alk  “2k  “21  “31  +  “lk  “21  “2k  “31  +  “U 


,  T|)  , 

^  e  *  Sk  e  6k+l 


.0^0, 


+  Sk(Sk  +  N2k^  e 


where  6^  is  the  Kronecka  delta  defined  by  5k_^  =  1  if  k 
if  k  f  2,. 


(E- 


(E- 


—  j(wkTj+V^) 

;32.  e 


“32,  “2k  “22, 


2,;  and  6k_A 


~tr*. 
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H  K  l  ,l-».  I  ** 


M 


-:v. 


■ 


► 


h: 


s 


k: 


m 


n 


Now  substituting  Equation  (E-ll)  into  Equation  (E-10)  yields; 


B  B 


ZI  (jt^)  (ju2)  hljc  h2j(  h?2  h32  Sj 


k=-B  l=~B 


u 

+  E  wk  hlk  h3k  lh2kl2  Sk 


k=-B 


-E 

k=-B 


“k  hlk  h3k  lh2kl2  VSk  +  N2k^ 


B 


B 


H  hlk  h2k  Sk||  L  H  h2k  h3k  Sk 
,k=-B  /  \  k=-B 


D 

-E 

k=-B 


fo  fo 
rl  r2 


“k  hlk  h3k  lh2kl2  Sk  N2k 


“k  hlk  h3k  lh2kt2  Sk  N2k 


( E-12) 


k=-B 


Finally,  using  Equations  (E-9a)  through  (E-12)  in  Equation  (E-8)  yields 


C0V(ti,t2)  -  —z 

I 

k*-B 


8 

k=-B 


w 

‘E  wk  hlk  h3k  lh2kl2  Sk  N2k 


8 

E  wk  hlk  h2k  Sk)f  E  “,2  ho'-  hr"'  S,‘ 

ik=-B 


( E-13 ) 


k  "2k  "3k  Jk 


E-7 
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Equation  (E-13)  can  also  be  expressed  in  continuous  frequency  domain  as 


C0V(Ti,t2)  =  -  ^ 


d 

J  to2  h^(to)  h^w)  | (a>)  j ^  S (to)  ^(oj)  do) 

2tt  _ 

T  B  B 

•  /•  2  *  /*  2  * 

/  oj  h^  h2 ( tu)  S(uj)  dco  /  a)  h 2 ( cu)  h^w)  S(cj)  du 

•'n 


(E-14) 


Note  that  if  either  the  definition  of  or  t ^  is  reversed,  the  sign  in 
Equation  (E-13)  and  (E-14)  is  also  reversed. 

Finally,  it  can  also  be  shown  similarly  that  the  variances  of  the 
estimates  are  given  by: 
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be  the  signal  and  the  noise  power  spectra  and  assume  identical  frequency 
spectral  filters  for  all  channels,  then  the  covariance  expression  of 
Equation  (E-14)  reduces  to 


E-9 
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C0V(^,t2)  = 


-1 


2 

7 


(E-17) 


and  the  variances  reduce  to 


VAR( )  =  VAR(t2) 


(S/n/ 

1  +Ts/n 


2 

Cl) 


2tt 

2  (S/N)2 

T 

3(1  +2  S/N) 

(E-18) 


Comparison  between  Equations  (E-18)  and  (E-17)  shows  that 


cov(^2) 


S/N 

•  r~rrm 


VAR(t,)  . 


( E-19) 


E-10 


APPENDIX  F 

DERIVATION  OF  EQUATION  (3.5.2-7a) 

From  Equation  (3.5.2-7a)  we  have 
B 

fr  ju)k  lhk  I2  Vk  fr  Vk  \  ’ 

k*l 

Using  the  relation 


8<P  _  3 $  5t1  3<S  9t2 

3r  3t^  3r  3t£  3r 

and  Equation  (3.5.1-10a-c) 


fe-  (vk  xh  vk*>  •  n  vk  % 


in  Equation  (F-l)  yields 


F-l 
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B 

r 

k=l 


lhkl2  ^  vk  vk 


( F-4) 


But 


Vk  Vk 


e  K  1  0  0 


0 

0 


1  0 
*Ht2 


0  e 


1  1  1 

JXti 

e  K  1  0  0 

1  1  1 

0  1  0 

1  1  1 

jXt2 

0  0  e  K  £ 

1 

%Ti 


-Hti  -n(i 

e  e 


"j(0kT2 


■H'W  J'V2 

e  e 


and 


«k*Vk  1„  Vf  -  |  l«lk  I  +  I 


=  !  k  j2  +  iot2k f2  +  |a'5'-12  ■ 


(  jtA-T, 

+  (a2k  alk  e  }  + 


3k 1  j 
(a2k  alk 
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Jr  *  ejV2w  ,  *  eJV2.*l 
j^a3k  a2k  ^  (a3k  a2k  ^  J 


j(a3k  alk 


eJWT2\  eJVTl+T2>  J 


)  +  (at,  a 


3k  lk 


)*  • 


( F-6) 


Now  using  the  fact  that  the  first  term  in  Equation  (F-6)  is  not  parameter- 
dependent  and  that  quantities  in  the  parentheses  are  conjugate  symmetric. 
Equation  (F-4)  can  be  rewritten  as 


3A(r,e) 


g 

If  2^  lhkl  (“2k  “lk  e  +  “3k  “2k  e  +  “3k  “lk  e 


oo 

Mw>|2  G(u;  Tl>  t2^  dw 


(F-7) 


where 


,  T  \t  *  \  ajVl  +  ,  *  ,  Ht2  +  ,  *  *  . 

3;  TV  t2)  ■  T  [(o^  alk)  e  +  a^)  e  +  («3k  «U)  e  J. 


( F-8) 


F-3/F-4 
Reverse  Blank 


APPENDIX  G 

CALCULATION  OF  LOCALIZATION  UNCERTAINTY  FROM  THREE-SENSOR  ARRAYS 

This  appendix  calculates  the  localization  uncertainty  based  on  a  two 
inter-array  time  delay  measurement. 

From  Equations  (D-lOa)  and  D-lOb),  the  approximate  (r  »  L)  time  delays 
to  range  and  bearing  relation  are  given  by 


—1  c 

e  =  -sin  {^j-  (ij  +  x2)} 


c 


2 

cos  e 
(t2  “t1> 


(G-l) 


(G-2) 


where  c  is  the  speed  of  sound,  and  L  is  the  inter-array  separation. 

T 

Letting  X  =  (e>  r)  ,  X  =  (Ti»  T2^»  Equations  (G-l)  and  (G-2)  can  be 
written  as 

Y  =  f (X) 

where  f(X)  =  [f^X)  1S  vector  function  defined  according  to  E- 

quations  (G-l)  and  (G-2).  Let  X*  and  X*  be  the  nominal  values,  6X  and  5X  be 
the  deviations  from  the  nominal;  then  a  first-order  expansion  about  the 
nominal  (R,  B)  yields 

3f* 

<S X  =  6X  ( G-3) 


where  it  is  defined  that 
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3f* 

3X- 


!2 

3f! 

3e 

3r 

3fJ 

3f| 

3e~ 

w~ 

(G-5a) 


where 


2L  cos  B 


( G-5b) 


b  = 


R2c 


L2  cos2  B 


(G-5c) 


since  R  »  L. 


Now  post-mu ltip ling  Equation  (G-4)  by  its  transpose  and  taking  the 
expected  value  yields  the  desired  covariance  relations: 


Let 


COV(X)  = 


3f*\ 

/ 

3X~  / 

C0V(X)( 

r 

o2 

o2 

11 

12 

o2 

a2 

12 


11 


(6-6) 


(6-7) 


then 


G-2 
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H 


COV(V)  = 


a2 

11 


a  2 
12 


a  2 
12 


a  2 
11 


-b 


„  2  .  2  v  2 

2(cth  +  o12)a 


2  2  -  2 
2(an  ‘  a12'b 


( G-8) 


and  the  determinant  of  COV(Y)  is 


|C0V(Y)|=  4  a2  b2(a^  -  a}2) 


(6-9) 


The  area  of  the  one-sigma  error  ellipse  is 
S  =  ir  |  C0V(  Y_)  1^ 


=  -2  ab  it 


V 


4 

CT11 


4 

a12 


(G-10) 


The  time  delay  covariance  matrix  for  the  optimum  approach  is  (see  E- 
quation  (3.5.2-16d)) 


VAR(^) 


(G-ll) 


G-3 
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and  the  conventional  is 


1 

S/N 

"  1  +  2  S/N 


S/N 

•  r  +  2  s/n 

1 


VARc(f1) 


(G-12) 


Thus  the  corresponding  area  is 


Tj.(.Rc)-_  vAR(t) 
LJ  cos'3  B 


( G-13) 


and 


TT  (  RC  )  ^ 

V  3 — 

LJ  cos  J  B 


V(1  +  3  S/N)(l  +  S/N)  VA  , 
(1  +  2  S/N)  VAVV 


( G-14) 


where  VAR(t^)  and  VARc(t^)  are  the  estimated  time  delay  variances  of 
the  optimum  and  the  conventional  processors,  respectively. 

Finally,  from  Equations  (G-8)  and  (G-ll),  the  optimum  range  and  bear¬ 
ing  variance  expressions  are  given  by 


VAR(B)  =  2(0^  +  o\2)  a2 


1 

”  7 


VAR(tx) 


( G-15 ) 


TR  6757 


i? 


Kf 


» 


KV 


* 


h 


i 


VAR(R)  =  2(a^  -  <^2)  b2 

’3c2(nJrTr)SAR<;i> 


and  from  Equations  (G-8)  and  (G-12),  the  conventional  range  and  bearing 
variance  expressions  are: 


( G-16 ) 


VARc(s!  ■  i  (rcfnrf  (lVrf/x) 

VflRc<R)  *  2c2(r4nrf  ( 1  *  1  !/n ) 


Using  the  relation  given  in  Equation  (3.5.2-16c), 


VARfrj)  _  2  /  1  +  3  s/N  \ 

'  i  (rrrw). 


the  conventional  range  and  bearing  variance  equation  can  also  be  expressed 
as: 


(G-17a) 


(G-17b) 


varc(b)  *  |  (rim) 2  »ar<V 

■  3  (tVt57h)  VAR(S) 


(6-18) 


and 


G-5 
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VARC(R) 


R 

L  cos  B 


VAR(^) 


=  VAR(R)  . 


(G-19) 


Equations  (6-18)  and  (G-19)  show  that  for  estimating  range,  the  conven¬ 
tional  approach  and  the  optimal  approach  have  identical  variance.  On  the 
other  hand  for  estimating  bearing,  the  conventional  yields  a  variance 
which  is  three  times  the  optimum  at  low  SNR  environment. 


APPENDIX  H 

GCC  STATISTICAL  PERFORMANCE  IN  THE  PRESENCE  OF  INTERFERENCE 


This  appendix  develops  the  appropriate  expressions  for  the  bias  and 
variance  of  the  time  delay  estimate  from  a  Generalized  Cross-Correlator  (GCC) 
in  the  presence  of  interference. 

Let  the  frequency  domain  representation  of  a  T-second  observation  of 
waveforms  from  two-sensor  arrays,  A^  and  A^,  in  the  presence  of  J  targets  be 
written  as: 


aik =  hk  h  +  nik 


a2k  "  -2k  £k  +  n2k  *  k=1’  2»  **•»  8 


(H-lb) 


where  B  is  the  highest  frequency  components  of  the  signal  or  the  noise  pro¬ 
cesses;  and  n2^  are  the  Fourier  components  of  the  noise  processes  at 

frequency  =  2Trk/T.  In  addition,  the  complex  vectors  B^,  V^k  and  are 

defined  by 


ilk  -  » 


J“k°il  H°12 
,  e 


HDiJ 

e  ^  1Jl;  1  •  1,  2 


(H-lc) 


and 


=  f6kl’  8k2’ 


( H-ld) 


where  D-.  is  the  propagation  time  delay  from  target  j  to  sensor  i  and  B.  .  is 

1 J  K  J 

the  Fourier  conponent  of  the  signal  spectrum  of  target  j.  For  convenience,  it 
is  assumed  that  both  the  signal  and'  noise  processes  are  zero  mean  and  mutually 
uncorrelated.  Thus,  we  have  the  relations  for  the  discrete  signal  and  noise 
power  spectra 


6ki  6kj 


skj;  1f  1  s  J 


0  ;  if  i  t  j 


( H— 2a) 


^ik  njk 


if  i  =  j 


0  ;  if  i  t  j 


(H-2b) 


where  (  )  denotes  the  expected  value. 

Without  loss  of  generality,  let  target  j  =  1  be  the  target  of  interest, 
then  the  best  estimate  of  =  D22  -  from  a  GCC  (see  Equation  (3.5.1 -14b ) ) 
is  to  seek  the  null  of  the  equation: 


f(T>  *  £  HIM2  °ik  “Ik =0 


(H-3a) 


2 

where  h^  ,  given  by 


,2  .  Skl/Nk 


(H-3b) 


is  the  optimum  spectral  shaping  filter  for  a  single  target  environment. 
Because  of  the  observation  noise.  Equation  (H-3)  is  a  stochastic  algebraic 
equation.  On  average,  however,  the  mean  of  the  estimate  must  satisfy  the 
equation 


7(t)  =  X  jwk!\l2  * 


lk  “Zk 


e  *  »  0. 


( H-4a) 
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But  from  Equations  (H-la)  and  (H-lb),  one  obtains 


alk  a2k  =  -ik  Sk  -2k 


j3l 


■no¬ 


where 


TJ  ■  °2j  -  Dlj 


is  the  actual  time  delay  difference  between  sensor  array  A 
and 


sk 


diag{Skl,  Sk2,  .  .  SkJ} 

is  a  diagonal  signal  spectral  matrix.  Therefore,  the  mean 
satisfy  the  equation 


ttft  -  jn  Hihki2  ^lk  sk  %  eJV=  ° 

k=-B 


or  using  Equation  (H-4b),  we  have 


J  B 


TRT  =2^2^  ja)klhkl  Skj  e 


J“k(T-Tj) 


j=l  k~-B 


Equation  (H-5b)  can  be  manipulated  to  yield 


( H-4b) 


(H-4c) 

and  A2  of  target  j 


(H-4d) 

estimate  must 


(H-5a) 


( H-5b) 


H-3 
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^  *  St  E  V  t 

j=l 


Tj] 


=  0 


where 


Vt)  'E  lhkl2SkJe 

k=-B 


jiat 


is  the  auto-correlation  of  target  j. 


Denoting  the  total  power  for  target  j  by 


u 

’j  ■  E  v 


k=-B 


the  total  received  power  by  the  sensor  array  as 
J  B 

p  -  2  pj  +  N;  N  =  E  nk- 

j=l  k=-B 

and  the  normalized  total  correlation  by 


p(t)  5^, 


then  Equation  (H-5c)  can  be  rewritten  as 


p(t)  -  2^  a.  ^Pj(T  -  Tj)  -  0 


j=l 


(H-5c) 


(H-5d) 


(H-6a) 


(H-6b) 


(H-6c) 


(H-7a) 
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where 


*  Yp 


(H-7b) 


and 


pj  (  t)  =  R j ( t) /Pj •  ( H-7c) 

In  general.  Equation  (H-7a)  must  be  solved  numerically;  however,  for  the  case 
in  which  the  interference  is  close  to  the  true  target  of  interest,  then  one 
can  write 


P'\(0) 

pj(°)  +  Pj(0)(T  -  Tj )  +  —  (T  -  Tj ) 


(H-8a) 


where  p'.(0)  and  p!?(0)  denote  the  first  and  second  derivative. 

J  J 


But 


juk|hk|2  skoVj  *  ° 

\k=-B  / 


(H-8b) 


and 


Pj(°) 


2  SkjJ/Pj 

>k=-B 


/  8 
-Alf 


w£|hwj^  S.(u)  du) 


U 


TR  6757 


T(Sj/N)(S1/N) 
iPj(l  +  2  Sj/NJ 


co2  dto 


(H-8d) 


where  in  the  last  equation  the  signals  and  noise  are  assumed  flat  and 
band  limited. 


Therefore,  using  Equation  (H-8a)  in  Equation  (H-7a)  yields 


y;  3j  p  j  ( o )  (  t  -  Tj)  =  o . 
j=i 

Thus,  the  mean  estimate  of  is 


\j=l  /  \j-l 

■(£"■)  (£'v 


(H-8e) 


where 


D 

w  j~\f  CD2 1  h  (a>)  1 2  Sj{co)  du 


(H-8f ) 


Now  using  Equations  (H-8e)  and  (H-7b)  in  Equation  (H-8d)  yields 


-  T7->  i  r j  i-i 

T1  ’  2^(S,/N) 

j=l 


(H-9) 


H-6 
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Thus,  the  time  delay  estimate  from  a  second  expansion  of  the  correlation  func^ 
tion  is  a  weighted  linear  combi  nation  of  time  delays  from  every  target.  For 
the  two-target  case,  let  the  signal  of  interest  be  =  S,  interference 
S2  =  I,  then  Equation  (H-9)  yields 


:  _  s  A  1 

T1  “  TT 7  T1  TTT  t2  • 


(H-10) 


The  bias  in  general  is 


bl  =  T1  -  T1 


^  (Sj/N)(Ti  -  TX)  {Sj/N) 

.j-i  J  L  j=i 


(H-ll) 


and  for  the  two-target  case 


h  =  r~  S/I  <T2  '  Tl)> 


(H  -12) 


We  next  derive  the  expression  for  the  variance.  Using  the  linearization 
procedure,  it  is  easy  to  verify  that  the  resulting  variance  of  the  time  delay 
estimate  from  seeking  the  null  of  the  function  f ( t)  is  given  by 


VAR<V  * 

(  8f<Tl>  ) 


(H-13a) 


where  Tj  is  chosen  such  that 


H-7 
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f(Ti> 1  2  niV2  “TT“2k eJVl  *  °- 


But  from  Equation  (H-3a)  and  (H-5a-d),  we  have 


x  B  B 

f(Tl)  =  S  Z  (j“k)( M)  lhkl2|h4I2  aik  “2k  au 

k=-B  £=-B 


3f(T,) 


~  =  £  (juk)2|hk|2  alk  ^k  e 


nTi 


Now  using  the  following  relations 


alk  “2k  =  -lk  h  ^2k  =  Sk.i  e  ^  J 


alk  °2k  aU  =  alk  “Sk  aU  a2A  +  alk  aU  ®?k  a2Jl  +  alk  “Jt  a2k  al£ 


alk  a2k  aU  “Ji  +  (-lk  Sk  -fk  +  V  6k+H 


+  dik  sk  \-i 


(H-15b) 


Equations  (H-14a)  and  (H-14b)  can  be  simplified  to 
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f(T, ) 


■  E  “kiM4k  K  iik +  v2  -  (iik  \  %  eJVl>2 1 


1;  k=-B 


(H-16a) 


37TrTT  2 

r-*-E  "k'h 

k=-B 


3t 


,2  u  :  JwkTl 

k  I  -lk  Sk  y$k  e 


(H-16b) 


Substituting  Equations  (H-16a)  and  (H-16b)  in  (H-13a)  yields  the  desired 

A 

expression  for  the  variance  of  the  Tj  estimate: 


VAR(Tl)  = 


B 

£ 

k=-B 


E  “k2ihki4!<iik  \  it 


2  ~  JOJT,  ?  ) 

(!ktNk>  -<iiksk%e  >1 


B 

l 

k=-B 


E  “lihki2(iu  5k  iik> eJVl 


1  2 


(H-17a) 


/ 


_  2tt  B 


2  A  -  2  ,  -  j^T-i  ,  2 

w  |hf  [(V1  S  VJ  +  Nr  -  (Vj  S  V|  e  1) 


do) 


(^*oj2|h|2  lx  S  V|  eJU5Tldo.) 


(H-17b) 


Recall  that 


V l  S  V|  +  N  =  N  + 


Esj 


(H-18a) 


V,  SV*e 


Jwt1 


■  y.s. 


jcoA. 
1  J 


( H-18b) 
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For  the  special  case  where 


S j ;  <_a)_<  1^2  and  -u^  _<  oa  _< 


0  ;  otherwise 


IN;  _<  a)  <  ^2  and  -wg  £  u  _<  -o>^ 
0;  otherwise 
Equation  (H-19)  reduces  to 


sV  R(0)  -EE  s.  s.  R(A.  +  A,) 
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k 


where 


ai  -  si 


( H-20b) 


co  cos  cox  dco 


( H-20c) 


p(t)  =  R(t)/R(0) 


(H-20d) 


1  ;  if  t  =  0 


(H-20e) 


r  2  ~\ 

1  2co  cos  cot  ,  (cot)  -  2  j  « 

7T7TT  - 5 -  +  - 1 — 0 -  Sin  COT  ,  IT  T  f  0. 


Woj 


For  the  two-tahget  case.  Equation  (H-20)  becomes 


*  1  -(a,  p(2Ai )  +  dp  P ( 2 A p)  +  2a, a?  p(A,  +  A«)]  , 

VAR(t,  )  - - - - - - - - - =-=-* — - - —  R-1  ( 0) 

2[ax  p(Aj)  +  a2  plfyr 


(H-21a) 


Note  that 


(H-21b) 


(H-21c) 
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where  we  have  let  *  S,  S2  =  I  be  the  signal  and  the  interference  power 
spectrum,  respectively.  If  the  interference  is  far  away  from  the  target  in 
time  delay.  Equation  (H-21)  becomes 


VARJ-^)  3  1im  VAR(x1) 


Therefore,  the  normalized  variance  is 


(H-22) 


VAR(tj) 

var„( ?i ) 


1  -  (a^  p(2Aj)  +  a^  p(2A2)  +  2a^2  p(A^  +  ^ 
2Uj  p(A:)  +  a2  P(A2)) 
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Appendix  I 

CALCULATION  OF  MULTITARGET  GENERALIZED  CROSS-CORRELATION  COVARIANCE 


The  calculation  of  the  covariance  matrix  resulting  from  discrete  obser¬ 
vation  of  the  Generalized  Cross-Correlation  (GCC)  function  in  the  presence  of 
multitarget  interference  is  presented.  The  results  presented  here  are 
obtained  from  a  frequency  domain  formulation. 

Let  the  noisy  observed  waveforms  from  two  sensors  in  the  presence  of  J 
targets  be  represented  by 

J 

yi(t)  =  23  s^t)  +  ni(t)  (1-1) 

i=l 


J 

y2(t)  Si(t  +  V  +  n2^  (1-2) 

i=l 

where  s^(t)  is  the  ith  target  signal  waveform  and  n^(t),  n2(t)  are  the  noise 
processes.  It  is  assumed  that  signals  and  noises  are  zero  mean,  mutually 
uncorrelated,  band-limited  Gaussian  processes. 

Let  the  waveforms  be  sampled  at  a  sampling  rate  of  At  seconds  such  that 
T  *  N  At  seconds  of  the  waveforms  are  observed.  Discrete  Fourier  transforms 
of  Equations  ( I -1 )  and  (1-2)  yield  the  equivalent  frequency  domain  representa¬ 
tion  as 

J 

°lk  ■  2  8ki  +  "Ik  d-3) 

i=l 


1-1 
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^k 


J 


i=l 


+  ^2k 


where 


=  2irk/T 


aAk  =  lSyA(nAt)  e"j(2imk/N)  ;  i  =  1,  2 
n=l 


6U  e-J(2mk/N) 


i  =  1,  2,  ...  c 


n=l 


and 


N 

"Mt  ’  S  S  "t(n4t) 


n=l 


e-j(2imk/N)  . 


I  =  1,  2  . 


The  GCC  is  obtained  from  the  following: 

R(t)  =  £  “lk  “Ik  lHkl2  ^ 

k=-B 

2 

where  |Hk|  is  the  spectral  shaping  filter. 

Equation  (1-8)  is  usually  implemented  via  inverse  FFT. 
discrete  observations  of  the  GCC  R(t),  R(nAt),  are  obtained, 
let  R(n)  =  R(nAt),  then  we  can  write 


(1-4) 


(I-J 


(I-< 


(1-5 


(1-8 


Consequently, 

For  simplicity. 
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B 


*<">  =  2  °ik 

k»-B 


jHk|2  eJ(2lTnk/N)  . 


(1-9) 


We  are  interested  in  the  auto  covariance  of  R(n)  and  R(m),  Anm*  for  any  n 
and  m.  In  general,  one  can  write 

Anm  =  RlnT^TmT  -  ^TTf  RM  •  (1-10) 

Now  substituting  Equation  (1-9)  into  ( I -10 )  and  simplifying  yields 


k=-B 


alk  “Ik  aU  *21  '  alk  “2k  aU  °?£  ) 


I  |2  e j ( 2tt/N ) (kn+fim) 

St 


I 

1 


(1-11) 


Since  a.^  for  all  i  and  k  are  zero  mean  complex  Gaussian  random  variables 
uncorrelated  for  different  frequency  one  can  write  (using  the  fourth  order 
product  moment  formula) 


“lk  ^k  °U  =  alk  “Sk  aU  +  alk  °l«,  “2 k  *$1  + 


“lk  °2i  °2k  ®U  * 


(1-12) 


Substituting  Equation  (1-12)  into  (1-11)  yields 


B  B 


nm 


s  ]C  2£{(°ik  an  °ik~ah  +  a!k  a2k  an) 


k=-B  f-B 


|Hk|2  IHjI2  ej(2^/N)(kn+tm)j_ 


(1-13) 


1-3 
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1-4 
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^2  (S2  +  2SkNk  +  N2)  |Hk|4  eJ'(27rk/N)(n_m) 


+  |Hk|4  S2  ej(2rrk/N)(n-Hn-2D/At) 


(1-17) 


□ 

~h  f  +  2SN  +  lH!4  eJ'^n'm^At  da) 

J  D 


*f  |H|4  S2  eJ"((n<”>“-20)  do. 

J  D 


d-18) 


where  S,  N,  and  H  denote  the  continuous  power  spectra. 


For  the  case  of  uniform  spectral  shaping,  i.e.,  |Hp  =  1,  we  obtain 
=  P<j$(^  -  w)  +  2p^(n  -  m)  +  p^(n  -  m)  +  P^^[(p  +  m)  -  2D/At] 


(1-19) 


where  if  $xx(u))  is  the  power  spectrum  of  X,  then  pxx(t)  is  given  by 


P*y(t)  =  Yfi  f  ^yy(uj)  eJaJT  dw  • 

Q 


(1-20) 


For  a  two-target  case  with  equal  noise  power  spectra.  Equation  (1-16) 
reduces  to: 


1-5 
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$  ’  E  <skl  *  Sk2  *  Nk»2  lHkl4  eJ(2"k/NH'’-m> 


t  (hi  ^ 


IHJ"  e 


4  j(2-rrk/N)(n+m) 


(1-21) 


B 

=  2  fHkl4  )  [Skl  +  Sk2  +  Nk  +  2(SklSk2  +  SklNk  +  Sk2V) 

i.  _  n  * 


J(2irk/N)(n-m)  +  $2  *j2ajkDl  2$  $  -H(D1+D2) 

[\l  e  "klTc2  e 


+  Sk2e 


‘j2a)kD2  e  j  (2irk  /N )  ( n+m) 


(1-22) 


Hence  the  continuous  approximation  is  given  by 


B  ( 

"fey*  lH!4  |  S2  +  S2  +  N2  +  2( SXS2  +  SjN  +  S2N)|  ejh,(n'm)At 

T  «  -j2o)Di  -jw(  Di+Do)  _  —  j 2a)Do 

+  [S2  e  ♦  ZSfy  e  ♦  S2  e  2, 


>(n+m)At  du  p 


(1-23) 


1-6 


2 

Now  for  uniform  spectral  shading  |H|  =  1  and  using  the  definition  of 

Equation  (1-20)  in  (1-23)  yields 

(2) 

=  ‘  m)  +  pS2$2{n  “  m)  +  pNN(n  '  m) 

+  2  ‘  m)  +  pS2N(n  '  m)  +  pS2N(n  "  m) 

+  PS1S2  +  M)  -  20-^/At ]  +  2Psis2  f^n  +  M)  "  (Dj  +  D 2 ) / At ] 

+  PS  s  f(n  +  M)  -  2D2/At]  .  (1-24 
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